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EVALUATION 


This  effort,  in  support  of  TPO  IC2a,  Aircraft  Identifica¬ 
tion,  was  a  study  to  assess  the  relative  effectiveness  of 
various  glint-error  suppression  techniques  in  radar  systems. 

Methods  of  pulse-to-pulse  decorrelation  investigated  in  this 
effort  included  frequency  diversity,  spatial  diversity, 
temporal  diversity,  and  polarization  diversity.  The  feasibility 
of  adaptive  processing  techniques  was  also  assessed.  This  study 
has  advanced  our  understanding  of  the  phenomenology  of  both 
range  and  angle  tracking  errors.  It  has  identified  the  significance 
of  both  types  of  errors  as  a  function  of  target  dimensions, 
geometry,  and  range.  Finally,  this  study  has  identified  the 
optimum  decorrelation  techniques  and  associated  signal  process¬ 
ing  techniques  for  minimizing  target  glint-induced  errors. 
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1.  INTRODUCTION  AND  SUMMARY 

An  analysis  of  target  noise  effects,  with  particular  emphasis  on 
target-induced  radar  positional  measurement  errors,  was  conducted.  Several 
new  and  significant  findings  were  derived  for  Rayleigh  targets  and  showr»  to 
be  applicable  to  complex  targets,  such  as  aircraft,  via  a  detailed  simulation 
analysis. 


For  such  targets,  the  target  noise  effects  were  shown  to  be  con- 

2 

trolled  by  the  target's  RCS  centroid  and  second  central  moment  ,  in  the 
radar  measurement  dimension  of  interest.  The  RCS  centroid  was  shown  to 
determine  the  mean  target  location  under  hi^  signal-to-noise  conditions. 
Assuming  N  independent  radar  measurements  are  averaged,  the  Cramer-Rao  bound 
on  the  ref'ultant  positional  variance  was  shown  to  be  W^/SN.  In  situations 
\'rfiere  it  is  practical  to  implement  a  weighted  average,  RCS  weighting  was 
sho',va  to  be  optimum  and  results  in  a  measurement  variance  of  /P?(N-1 ) .  For 
those  applications  in  which  such  weighting  is  impractical,  data  editing  with 
an  RCS  threshold  of  six  tenths  of  the  mean  target  RCS  was  shown  extremely 
effective  and  results  in  a  variance  of  1.5wf/2N. 

For  non-Rayleij^  targets,  such  as  a  simple  two  scatterer  target, 

RCS  weighting  was  shown  to  be  effective  but  not  optimuin.  While  this  was  not 
pursued  to  any  great  extent,  a  heavier  RCS  weighting  appeared  preferable. 

In  such  situations  a  technique  previously  reported  in  the  literature*  proved 
better  than  RCS  weighting. 

This  report  is  organized  into  two  principal  sections.  Section  2 

addresses  the  impact  of  target  effects  on  the  noncoherent  detection  problem 

associated  with  search  radars.  Generally  speaking,  target  detTOtion  was 

shown  to  degrade  with  increasing  pulse-to-pulse  correlation.  Temporal  and 

frequency  diversity  techniques  for  decorrelating  the  target  signal  were  sliown 

2 

to  depend  on  the  target  parameter  \  in  the  cross-range  and  range  dimension 

♦  Loomis,  J.M.  and  Graf,  E.R.,  "Frequency-Agile  Processing  to  Reduce 
Glint  Pointing  Error,"  Trans.  IEEE,  vol,  AES-10,  Nov.  1974. 

pp.  811-820. 
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respectively.  For  representative  operational  conditions,  frequency  diversity 
was  shown  to  be  the  preferable  technique. 


Section  3,  vAilch  deals  with  the  target  location  studies,  represents 
the  major  study  emphasis.  Analytic  derivations  and  computer  simulation 
validation  of  the  key  results  cited  earlier  are  presented. 
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2.  TARGET  DETECTION  STUDIES 


In  this  section  the  role  of  waveform  diversity  in  enhancing 
fluctuating  target  detection  is  addressed.  The  principal  finding  is  that 
frequency  diversity  rather  than  temporal  diversity  is  the  most  effective 
means  of  enhancing  target  detection. 

2.1  Target  Model 

2.1.1  Amplitude  Characterization.  Since  the  target  detection  study  was 
analytic  in  nature,  an  analytic  description  of  a  representative  target  radar 
cross  section  (RCS)  model  was  required.  In  this  section  a  very  brief  motiva¬ 
tion  and  description  of  the  selected  target  model  is  presented. 

The  targets  of  interest  to  this  study  are  principally  aircraft 
which  are  highly  complex  targets  having  many  scattering  centers.  For  such 
targets  a  Rayleigh  target  amplitude  fluctuation  model  is  usually  quite  repre¬ 
sentative  and  was  assiomed  for  this  study.  VAxile  this  distribution  results 
from  the  assumption  of  an  infinite  number  of  equal  scattering  centers,  tlie 
assumption  is  really  quite  representative  when  the  number  of  near  equal 
scattering  centers  is  as  small  as  4  to  6. 

This  result  can  be  justified  by  observing  the  data  presented  in 
Figure  1.  Here  the  target  RCS  was  assumed  to  be  the  sum  of  N  independent 
equal  amplitude  scatterers.  normalized  such  that  the  average  RCS  is  unity,  i.e. 

"""n  ■  I  ‘ 

where  are  independent  random  variables  uniformly  distributed  over  (0,271  ). 
The  cumulative  distribution  plotted  in  the  figure  is  * 

P(x)  =  P(RCS|^*x) 

=  y)dy 

*  Greenwood,  J.A.,  "The  Distribution  of  Length  and  Components  of  the  Sum 
of  n  Random  Unit  Vectors",  Annals  of  Mathematical  Statistics,  Vol.  26, 

No.  2,  June  1955. 


3 

I 


666* 

6666 


- 


Here  one  notes  that  fr^r  N  equal  to  4  to  6,  a  good  fit  to  the  Rayleigh  distri¬ 
bution  (N  =  00  )  is  achieved.  Since  it  seemed  plausible  that  a  generic  A/C 
model  would  contain  at  least  4  to  6  comparable  scattering  centers,  the 
Rayleigii  target  model  was  selected.  (Note:  In  Section  3.2.1  the  cross  sec¬ 
tion  distribution  resulting  from  a  detailed  aircraft  model  is  shown  to  be 
effectively  Rayleigh  for  either  horizontal  or  '-ertical  polarization). 

2.1.2  Spatial  Characterization.  Before  the  effects  of  diversity  on 
target  detection  can  be  addressed,  the  target  must  also  be  described  in  terms 
of  its  spatial  distribution  as  illustrated  in  Figure  2.  This  characterization 
will  also  be  necessary  for  the  theoretical  assessment  of  target  induced  radar 
measurement  errors  such  as  glint  (see  Section  3.2).  Here  the  complex  scatter¬ 
ing  amplitude  arising  from  a  pcint  (x,y)  on  the  target  of  Figure  2  is 


Figure  2  Rough  Radar  Target 


5 


represented  by  a(x,y  ).  Further  assume  that  scattering  f^om  any  differ¬ 
ential  element  is  totally  Independent  from  any  other  element,  has  random 
phase  and  the  expected  power  can  be  defined  by  a  function  s(x,y  ),  where: 


E  [  a(x3^,yj^)a*(x2,y2)  J  =  s(xj^,yj^)  5  (x2-Xi,y2-yi) 

Here  S  represents  the  Dirac  delta  -  function. 

If  the  target  extremes  are  denoted  by  +Wx/2  and  +«/y/2,  then  the 
complex  narrowband  return  signal  at  frequency  f  can  be  represented  as 

^y 

/1(f)  =  EXP  (-j47r  R^f/c)  j  j  a(x,y)EXP(-j4Tf  yf/c)dxdy 

-^x  -^y 

This  signal  is  assumed  to  be  a  circular  (i.e.  complex)  Gaussian  process  since 
it  is  the  sum  of  many  random  elements.  As  a  result  its  amplitude  distribution 
is  Raylei^,  vrtiich  is  in  keeping  with  the  assumptions  of  Section  2,1.1. 

With  this  signal  representation  the  correlation  parameters  of  the 
process  in  frequency,  time  or  space  can  be  addressed.  First  consider  tiie 
correlation  between  two  signals  at  different  centrr  frequencies.  Defining  the 
correlation  coefficient  as 

E  [/I(fi)//*(f2)] 

E[M(fi)l  2] 

J j  s(x,y)EXp[-j4Try(fj^-f2)/cjdxdy 
j j  s(x,y)dxdy 

and  a  normalized  range  scattering  function  Sy(y)  as 
Sy(y)  = 

yjs(x,y)dxdy 


the.  fi'equency  correlation  coefficient  becomes 
pH.::  •  y  S^(y)E:<P(-j2llWy)>:^y 

piO)  +  3p'(0)  + 

1  -  ^BjySy{y)dy  -2n^B^jy^S„iy}dy 

for  small  values  of  iS  =  2(f^  -  f2)/c.  Assuming'  the  target  measurements  were 
made  at  a  range  corresponding  to  the  range  RCS  centroid  of  the  target, 
then  yD'(o)  is  zero  and  p(B)  is  real  and  approximated  by 

p(B)  =  1  -  i>rr^/?^|y^Sy(y)c'y 

As  a  result  the  correlation  coefficient  in  frequency  is  seen  to  be  controlled 
by  the  second  central  moment  of  the  normalized  range  scattering  function. 

In  an  analogous  manner  the  correlation  coefficient  for  target 
i\>i.ar.ion  can  also  be  derived.  Here  two  target  samples  are  assumeo  taken  with 
the  second  delayed  in  time  by  an  amount  t.  During  this  time  interval  the 
target  is  assumed  to  rotate  through  an  angle  In  this  case  one  finds: 

pl^)  ^  1  -  X 

where  S  (x)  is  the  normalized  cross  range  scattering  function  and  J  =  2f(U  t/c. 

•X  K 

Here  the  angular  or  temporal  correlation  function  is  seen  to  be  controlled  by 
the  second  moment  of  the  normalized  cross  range  scattering  function.  Details 
of  this  derivation  are  contained  in  Appendix  A.l. 


2.2  Effects  of  Signal  Correlation.  In  this  subsection  the  effects  of 
signal  correlation  cn  target  detection  are  addressed.  A  Rayleigh  target 
model  with  varying  pulse-to-pulse  correlation  is  assumed.  The  teirget  statis¬ 
tics  then  lie  between  those  of  a  Swerling  I  model  for  unity  pulse-to-pulse 
correlation  and  a  Swerling  II  model  for  zero  pulse-to-pulse  correlation. 
Results  are  presented  for  a  scanning  search  radar  which  is  assumed  to 
noncoherently  integrate  N  pulses.  Generally  speaking,  detection  performance 
is  shown  to  be  a  maximum  when  the  pulse-to-pulse  correlation  is  a  minimum. 

As  a  result  noncoherent  search  radars  should  routinely  utilize  pulse-to-pulse 
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frequency  agility  to  decorrelate  target  statistics  and  thereby  maximize 
target  detection. 

For  a  Rayleigh  target  the  following  expressioriti  for  the  probability 
of  detection  and  false  alarm  were  derived  (see  Appendix  A. 2}: 

M 

■  JZ 

m=l 

PpA  =  EXPC2^f/N^) 

m=0  ro  J 

where  M  =  Number  of  Pulses  lntet’ra'>^£>d 

^  =  Threshold  Level 

Nq  =  Noise  Power 

/ifn  =  The  Eigenvalue  of  the  Received 

Signal  Covariance  Matrix 

The  expression  fcr  Pp  assumes  that  the  eigenvalues  are  distil*^- 1  and  must  ue 
modified  to  handle  the  case  of  repeated  eigenvalues. 

( 

Figure  3  is  a  plot  of  P^  vs  the  interpulse  correlation  coefficient 
for  M  =  2,4  and  o  pulses  and  signal-to-noise  ratios  of  3,6  and  9  db  wiieii  the 
false  alarm  probability  is  Aside  frc*^.  the  cast  of  M  =  2  with  3NR^3  db, 

detection  performance  is  seen  to  iuo'-cases  C.ocrvuon^g  a  -orp  ilse  corre¬ 

lation.  This  trend  is  increasingly  strong  with  increasing  signal- to-noise 
ratios. 

2.3  Operational  Implicat i ons .  Having  derived  expressions  for  tiie  pulse- 
to-pulse  target  correlation  in  Section  2.1.2  and  demonstrated  in  Section  2.2 
that  detection  performance  is  maximized  tdien  the  pulse-to-pulse  oorrelati  vi 
is  minimized,  we  address  the  operational  implications  of  these  findiiig.r.  :n 
this  section.  Frequency  agility  is  shown  to  be  the  most  pract  ical  ni  ums 
decorrelating  the  received  signal. 


Assuming  a  pulse-to-pulse  target  correlation  coefficient  of  .5  is 
desired,  the  expression  for  p(BI  derived  in  section  2.1.2  may  be  set  equal 
to  .5.  Noting  B  =2b£jc,  one  may  solve  for  A  f  as 


where  is  the  second  central  moment  of  the  target’s  RCS  function  in  range. 
This  relationship  is  graphed  in  Figure  4.  Indicated  on  this  figiore  are 
typical  values  of  Wp  in  meters  for  both  a  representative  cruise  missile  and 
fighter  aircraft.  These  values  were  derived  from  the  target  models  discussed 
in  Section  3.3.1.  For  these  values,  pulse-tc-pulse  frequency  agility  re¬ 
quirements  ranging  from  4.5  to  15  MHz  are  noted.  These  values  are  well  within 
the  capabilities  of  current  radars. 


To  determine  temporal  decorrelation  requirements,  the  temporal 
correlation  function  derived  in  Section  2.1.2  is  also  set  equal  to  .5,  This 
yields 


vdiere  A  is  the  radar  wavelength,C»)f{  is  the  angular  rotation  of  the  target 
relative  to  the  radar  line-of-sigjit,  and  is  the  second  central  moment  of 
the  teu7get’s  RCS  function  in  the  cross-range  dimension.  This  relationship  i:; 
shown  in  Figure  5  for  a  target  rotational  rate  of  1°  per  second  and  radar 
wavelengths  of  .01  and  .1  meters.  The  values  of  Wx  indicated  are  also  baaed 
on  the  target  models  discussed  in  Section  3.3.1.  The  resulting  temporal 
decorrelation  times  are  seen  to  exceed  ttie  interpulse  period  of  representative 
search  radars  which  typically  run  several  hundred  pulses  per  second.  As  a 
result  the  temporal  decorrelation  of  targets  will  usually  be  too  slow  to 
provide  uncorrelated  returns.  This  temporal  decorrelation  should  be  adequate 
to  provide  scan-to-scan  decorrelation. 


Assuming  a  pulse-to-pulse  target  correlation  coefficient  of  .5  is 
desired,  the  expression  for  derived  in  section  2.1.2  may  be  set  equal 

to  .5.  Noting  B  =2^f/ct  one  may  solve  for  A  f  as 


Af  = 


where  ig  the  second  central  moment  of  the  target's  RCS  function  in  range. 
This  relationship  is  graphed  in  Figure  4.  Indicated  on  this  figure  are 
typical  values  of  in  meters  for  both  a  representative  cruise  missile  and 
filter  aircraft.  These  values  were  derived  from  the  target  models  discussed 
in  Section  3.3.1.  For  these  values,  pulse-to-yulse  frequency  agility  re¬ 
quirements  ranging  from  4.5  to  15  MHz  are  noted.  These  values  are  v\rell  within 
the  capabilities  of  current  radars. 


To  determine  temporal  decorrelation  requirements,  the  temporal 
correlation  function  derived  in  Section  2.1.2  is  also  set  equal  to  ,5.  This 


yields 


where  A  is  the  radar  wave length, CJp  is  the  angular  rotation  of  the  target 
relative  to  the  radeir  line-of-sight,  and  is  the  second  central  moment  of 
the  target's  RCS  function  in  the  cross-range  dimension.  This  relationship  i£ 
shown  in  Figure  5  for  a  target  rotational  rate  of  1°  per  second  and  radar 
wavelengths  of  .01  and  .1  meters.  The  values  of  Wx  indicated  are  also  based 
on  the  target  models  discussed  :.n  Section  3.3.1.  The  resulting  temporal 
decorrelation  times  are  seen  to  exceed  the  interpulse  period  of  representative 
search  radars  which  typically  run  several  hundred  pulses  per  second.  As  a 
result  the  temporal  decorrelation  of  target:s  will  usually  he  too  slow  to 
provide  uncorrelated  returns.  This  temporal  decorrelation  should  be  adequate 
to  provide  scan-to-scan  decorrelation. 
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Lr^ure  4  Frequency  Diversity  Requirements 


3 


TARGET  LOCATION  STUDIES 


3,1  Location  Measurement  Error  Components.  Before  addressing  the  target 
noise  component  of  the  radar  measurement  error,  it  will  be  instructive  to 
first  review  all  the  components  of  the  location  measurement  error.  Since  tlie 
equations  of  motion  which  describe  the  dynamics  of  a  target  under  track 
usually  apply  to  its  center-of-gravity  (c.g. ),  it  is  useful  to  reference 
measurement  errors  to  the  target's  c.g.  If  x  is  the  measured  target  location 
(in  range  or  angle  and  x^g  is  the  location  of  the  target's  center  of  gravity, 
the  measurement  error  is  then  defined  as 


=  Xg  +  Xn  +  xt 

Here  the  error  is  portrayed  as  having  three  components.  The  first 
component,  Xg,  is  the  limiting  system  error,  or  in  other  words  the  jueasurement 
error  that  v;ould  occur  in  determining  the  location  of  a  constant  point  target 
with  infinite  SNR.  This  error  results  from  system  hai-dware  linitatious  an,; 
usually  modeled  as  a  zero  mean  error  whose  standard  aeviation  is  somewhere 
between  1/50^^  and  1/25*'^  of  the  basic  radar's  measvare.Ment  resolution,  H, 

The  second  term,xn,  represents  the  additional  measurement  error 
associated  with  a  point  target.  This  error  is  due  to  both  thermal  r.cist  (.•1.*  .  ' 
and  taiget  fxuctuation  effects.  As  witli  the  system  aolae  term,  ito  mean  vaiu  . 
is  also  Zero.  For  a  constant  target  the  CiOir.cr -Rao  bouj'.iJ  ..n  the  var'ianC'.  cl 
this  term  is 


Here  the  resolution  parameter  R  is  ^/D,  where  A  is  tiie  radar  wavelci.gth  and  H 
the  array  diameter,  for  angle  measurements,  or  C/2BW,  where  C  is  the  speed  of 
light  and  BW  the  pulse  bandwidth,  for  range  measurements. 

When  the  target  RCS  is  not  constant  but  fluctuates  witli  a  probal  ilitv 
density  function  p(RCS),  an  average  Creuner-Rao  bound  may  be  computed  .is 


5^  _  1.5  dRCS 

"  j'p(ncs)  dRCS 

T 

If  the  lower  intergration  limit  is  taken  as  zero  for  the  Raj^leigh  targot 
model  of  interest,  then  this  quantity  is  unbounded.  In  practice  however  one 
requires  a  detectable  signal  before  a  measurement  is  made.  As  a  result,  the 
lower  limit  can  be  set  equal  to  a  RCS  threshold  value  of  -NolnPp^^  where 
Is  the  design  false  alarm  probability.  When  this  is  done  one  obtains 


Pp^  )  /(1+SNR)J 


1.5  r"  eJ-  lri(  PpA  > 

ff^d+gNR)  ?d[“  ^Fa  ^ 


vdiere  is  the  exponential  integral  and  Pq  is  the  detection  probability,  i.e. 
Pp  =  EXP^lnPpy^/(l+SNR)j  .  The  resulting  normalized  measurement  standard 
den'iation  is  plotted  in  Figure  6  for  Pp^^  =  lO"'^  and  Also  plotted  for 

comparison  is  the  Cramer-Rao  bound  for  a  constant  point  target.  Here  it  is 
observed  that  the  measurement  accuracy  for  a  Rayleigh  target  is  a  factor  of 
two  poorer  at  a  30  db  SNR. 


The  remaining  term,  X  ,  represents  the  measurement  error  introduced 
by  a  non-point  target  and  is  the  principal  term  of  interest  to  this  study. 


3.2  Theoretical  Results.  In  this  section  an  analytic  treatment  of 
target  noise  effects  on  target  location  measurement  accuracy  is  presented. 

The  intent  of  this  analysis  was  to  develop  a  mathematical  understanding  of 
target  noise  effects  upon  which  candidate  target  noise  reduction  techniques 
could  be  derived.  While  the  analysis  presented  is  for  monopulse  errors, 
analogous  results  may  also  be  derived  for  range  measurement  errors.  In 
either  case  a  number  of  highly  useful  results  were  derived  including: 

•  Expressions  for  the  mean  and  variance  of  targ.et 
noise  measurement  errors, 

.  Definition  of  the  optimum  data  weigliting/editing 
for  track  data  smoothing  in  the  presence  of 
limiting  target  noise. 


NORMALIZED  STANDARD  DEVIATION 


These  results  sei^ved  as  a  basis  for  the  target  noise  reduction  techniques 
used  in  the  simulation  analysis. 


3.2.1  Preliminaries.  The  analysis  begins  with  the  definitions  of  the 
target's  angular  scattering  function,  0((6),  as  a  sample  function  of  a  r sro 
mean  circular  Gaussian  process  over  the  angular  width  of  the  target,  i.e. 

0((e)  =  0  for  |©1>J^L 

The  autocorrelation  function  of  this  process  is  defined  by 

Ejot(ei)ol*(©2)]  =  S(©i)S  -  ©2) 

i.e.,  the  process  is  vrtiite  in  angle  with  expected  power  density  in  angle  given 
by  S(©). 


Modeling  the  two-way  voltage  response  of  a  pair  of  squinted  ampli¬ 
tude  momoDulse  beams  as  ^  a  ^  b  ©,  the  target  returns  may  be  expressed  as 

=j«<.(©)(a  +  b©)  d©  +  n^ 

•V* 

Here  the  receiver  noise  contributions  n^  are  zero  mean,  uncorrelated  circular 
Gaussian  processes  with  spectral  density  Nq/2.  Since C((©)  and  n^  are  both 
Gaussian,  the  returns  r^  and  any  linear  combination  of  ^hem  are  also  Gaussian. 
When  the  angle  between  the  squinted  monopulse  beams  is  eq'Jal  to  the  3  db  beam 
width,  ©BWt  then  a  =  {?/2  and  b  =  1.17/©gy. 

Defining  the  monopulse  sum  and  difference  channels  as  sums  of 

Gaussians 

Z  =  (r+  +  r_)/2a 
A  =  (r+  -  r_)2/b 

\4ilch  must  also  be  gaussian  with  their  probability  density  function  given  as* : 

p(5,A)  =  (2Tr)“^lAH  EXP(-3^(  5*  A*)A"^(4)) 

*  H.L.  Van  Trees,  "Detection,  Estimation,  and  Modulation  Theory"  Part  III 
John  Wiley  &  Sons,  Inc,,  1971,  Page  600, 


\ 


1 


where 


and 


2  WMu 


JXi,  -  “ 

j  3(0) de  + 

p  -  e{5A*|  - 

J  OS(0)d0 

fe^  se)de 

No/2a‘ 


As  a  result  the  statistics  of  the  sum  and  difference  channel  are  completely 
determined  by  the  noise  level  and  the  first  three  moments  of  the  target's 
angular  RCS  function  S(S). 

Before  continuing  it  will  be  useful  to  define  two  key  parameters. 
The  first  is  the  quantity,  m,  defined  as 


m  = 


Ml 


At  high  signal-to-noise  ratios  this  quantity  approaches  the  target's  anguitir 
RCS  centroid  location.  The  second  p>arameter  of  interest  is  defined  by  W 
where 

‘  2 


= 


M I 


Ml 


This  quantity  is  seen  to  approach  the  2'^'^  central  moment  of  the  target's 
angular  RCS  function  at  high  SNR  and  hence  is  a  measure  of  the  angular  w^dtli 
of  the  target. 


With  these  definitions  and  knowledge  of  p(Z,a)  a  number  oi  useful 
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expressions  pertaining  to  monopulse  errors  can  be  derived.  Defining  the 
complex  monopulse  estimate  as 


Z 

where 


and 


X  +  j  Y 


X  =  Re  [z] 
Y  =  Im  [z] 


monopulse  signal  about  boresite 
quadrature  monopulse  signal  about  boresite 


U  =  X  -  m 


monopulse  signal  about  the  angular 
location  m 


one  finds  (see  Appendix  A. 3) 
p(U,Y)  = 


ir[u2  +  y2  +  w2]2 


p(U) 

p(Y) 


_ w2 

2  [u2  ♦  w2 

_ w2 

2  [  Y^  +  ^ 


These  probability  density  functions  indicate  that  the  statistics 
of  both  U  and  Y  depend  only  on  the  parameter  W  and  not  the  target  location. 
Further  it  is  shown  in  the  Appendix  that: 

E(U)  =  E(X)  -  m  =  0 

E(Y)  =  0 

In  other  words  the  expected  value  of  the  monopulse  estimate  x  is  the  location 

m. 


3.2.2  Variance  Expressions  For  Data  Editing  &  Smoothing.  Tech¬ 
niques  for  combining  several  target  position  measurements  in  a  manner  which 
minimizes  the  target  positional  variance  are  of  greet  interest.  Generally 
two  distinct  classes  of  techniques  can  be  defined.  The  first  are  data 
editing  schemes  vdiich  "edit  out"  or  reject  some  raw  measurements  prior  to 


track  smoothing.  These  techniques  are  suited  to  track  smoothing  algorithms 
which  treat  all  input  data  equally.  The  second  class  of  techniques  are  those 
which  are  designed  to  work  with  tracking  filters  which  permit  input  data 
weighting.  Naturally  these  techniques  should  be  capable  of  near  optimum 
performance. 


Consider  first  a  technique  for  data  editing.  As  formulated  in  the 
previous  section,  the  variances  of  both  the  real  and  quadrature  monopulse 
components  are  equal  and  unbounded.  Finite  values  for  the  variance  can  be 
obtained  by  treating  the  subset  of  target  measurements  associated  with  de¬ 
tectable  signal  levels.  Since  this  is  \«^at  any  radar  really  implements,  such 
a  constraint  is  realistic.  The  statistics  of  the  monopulse  estimate  U 
conditioned  on  a  received  signal  strength  of  tx  =  | Z I  is  Gaussian,  i.e. 


The  variance  of  U  conditioned  on  Ot  is  then 

0(/|«  =  /z'^  “ 

ot 

Here  it  is  seen  that  this  conditional  variance  is  proportional  to  the  square 
of  the  parameter  W  and  inversely  proportional  to  the  received  signal  strengt' 
cx.  ,  a  result  large  glint  variances  are  associated  with  small  signal 
levels.  Assuming  the  receiver  has  a  signal  detection  threshold  of 
can  compute  the  average  glint  variance  associated  with  detected  signals  as 


_ 

I  p(»-)dfli. 


\«here 


p(«*  )  =  ^  EXF(-« )  (This  follows  since  I?  I  is  Raylaigh) 


and  'm  is  the  exponential  integral.  Figure  V  depicts  the  dependence  oi 
this  expression  on  x  =  ^'/jULz-  Usually  the  threshold  is  set  to  yield 
a  given  false  alarm  rate,  i.e.  CX.j  =  -Nq  this  case  x  becomes 

equal  to  -  .in(Ppy^)/(l  +  SNR).  For  PpA  =  10"®  and  S1«R  =  25  db,  x  becomes  eq'aa 
to  .0436  and  is  found  equal  to  2,7. 


NORMALIZED  THRESHOLD  -  x 


I 


where  x  =  IjJif.  •  The  first  term  represents  the  average  Cramer-Rao  bound 
for  a  Rayleigh  point  target.  This  term  v/as  discussed  in  Section  3.1  and 
plotted  in  Figure  6.  The  second  term  in  the  above  variance  expression 
represents  the  impact  of  a  finite  target  extent.  The  parameter  Wq  is  the 
second  central  moment  of  the  target's  angular  RCS  function  and  hence  is  a 
measure  of  the  target  extent.  For  example  if  the  target  RCS  was  uniform 
across  a  target  extent  0^,  then  would  be  equal  to  0^,,/  j~12.  For  illus¬ 
trative  purposes,  Figure  S  shoves  the  dependence  of  this  combined  monopulse 
standard  deviation  vs  SNR  for  “  0,  .025,  .05,  .075,  .1  and  a  Pyf^  OF  10“®. 

For  comparison  the  Cramer-Rao  bound  for  a  constant  target  is  also  presented. 
This  is  the  value  which  is  usually  used  to  predict  radar  measurement  perform¬ 
ance. 


Now  assume  one  has  a  sequence  of  M  independent  monopulse  estimates, 
each  of  which  passed  a  detection  threshold.  One  standard  method  of  combining 
such  data  is  to  simply  form  an  unweighted  average,  i.e. 


The  resulting  variance  is  then 

Tiie  reference*  contains  a  reprint  of  an  IEEE  paper  which  concludes  on  the  bra- it 
of  simulation  that  a  superior  technique  is  to  select  only  that  monopulse 
estimate  associated  with  the  largest  signal  strength.  For  the  Rayleigh  model 
assumed  in  this  section  it  was  possible  to  derive  an  expre.ruion  for  ihis 
rofarenced  technique,  namely 

N-1 

- Jii -  EXP  [{ N-X  El  r ( N-k 

k!{N-l-kl!  *•  I 

k=0 


^SLA 


N-l-k 


These  two  methods  (ure  compared  in  Figure  0  for  a  Uircshold  value 

*  "Radars,  Vol.6,  Frequency  Agility  and  Diversity."  edited  by  David  K.  B-urton, 
Artech  Publishers,  Also  see:  J.M.  Loomis  and  E.R.  Graf,  "Frequency-Agile 
Processing  to  Reduce  Radar  Glint  Pointing  Error,"  IEEE  Transactions  on 
Aerospace  A-  Electronic  Systems,  Tlovenber  1974, 
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corresponding  to  =  10“®  and  SNR  =  25  db,  i.e.  **  ■-^(PFA)/SrJR=.043. 

I  VAien  the  number  of  measurements  N  is  between  2  and  4,  the  Select  Largest 

I 

Amplitude  (SLA)  technique  is  slightly  better  than  an  unweighted  mean  average. 
For  N>  6  however,  the  SLA  technique  becomes  increasingly  poorer  than  the  mean 

i 

average.  As  a  result  computing  a  mean  average  is  indeed  a  reasonable  technique. 

By  increasing  the  threshold  value  beyond  the  detection  thres’iol 
the  variance  made  smaller  as  indicated  in  Figi.u'e  7.  In  doir.j, 

so  however,  the  signal  detection  probability,  Pp  =  EXP("  ).  al.o 

reduced.  In  certain  situations  it  is  possible  to  define  near  optimum  values 
for  0<T  •  Consider  the  tracking  example  in  which  one  is  concerned  •.•.■itl  the 
variance  associated  with  the  smoothed  target  position,  Q"r  .  VJhen  the  tai-  ,e'; 
motion  can  be  described  by  a  polynomial  of  degree  K,  it  is  v/ell  know  that 
the  target  position  variance  approaches 

for  large  M 

1  for  K  =  0 
4  for  K  =  1 
9/4  for  K  =  2 

Total  number  of  independent 
measurements. 

Constant  measurement 
variance 


To  study  the  effects  of  threshold  setting  on  ^  ,  it  will  be  asuu.ico 


that  the  tracking  filter  being  utilized  does  not  permit  input  data  wei,,lvting, 
i.e.  all  measurements  are  treated  equally,  and  target  noise  effects  dominate 
the  measurement  error.  A  binary  weighting  will  be  permitted  in  whicii  the 


At  this  point  it  is  natural  to  ask  whether  thresholding  some 
variable  other  than  the  received  signal  strength  would  yield  a  better  binary 
weighting  for  track  veiriance  reduction.  In  an  attempt  to  answer  this  question, 
the  utility  of  thresholding  the  quadrature  monopulse  signal,  y  •  v.’as  also 
Investigated.  Here  the  conditional  monopulse  variance  v/as  i^ound  to  be 


= 


(See  Appendix  A. 3) 


This  result  states  that  large  monopulse  variances  are  associated  with  large 
values  of  Y ,  As  a  result,  discarding  measurements  when Y  exceeds  some 
threshold,  is  a  reasonable  technique  to  investigate.  VThen  this  v/as  done 
the  resulting  track  variance  v;as  found  to  be  (see  Appendix  A. 3) 

(^^(Y^)  =  - 1 2^1  +  (w/Yt)^]-Ai[yt/v;  +Vi  +  (Y3./w)2]'| 

The  term  in  |  ]  was  found  to  have  a  minimum  value  of  3.43  for  Y-p/W  =  1.3. 

For  this  sitiiation  75%  of  the  measurements  are  accepted  and  the  minima;.; 
variance  becomes 


MIN 


3.4 


Comparing  this  result  with  that  obtained  for  thresholdin 


strength,  i.e.  1,5  clearly  indicates  that  editinr 


is  superior. 


g  the  received  signal 
on  signal  strength 


Now  consider  the  class  of  smoothing  techniques  v/hich  permit  var-Vjhl' 
weighting  of  the  radar  measurements.  For  Rayleigh  target,  a  hirM  -"'-hd  sclut  i 
for  the  maximum  likelihood  estimate  of  the  target's  RCS  centroid  v/ar,  lour.*,  t  . 


be 


ML 


n=l 


rfcl 


(See  Appendix  A. 5) 
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1  .  poLjilion  estimates,  X^,  should  be  weighted  by  the 

assoc  ovLed  received  si-;nal  strength.,  0(j^  =  jZ  nP*  variance  of  this 

estimator  was  determined  (see  Appendix  A.b)  to  be 

2 

Comparing  this  result  with  those  for  the  unweighted  average  and  the  SLA 
technique  depicted  in  Figure  0,  the  PCS  weij^Jiting  is  seen  to  be  superior 
for  IJ>1.  (Mote  the  apnarej-tt  discrepency  at  N-1  arises  because  no  detection 
t’-resholdin,;,  was  assumed  in  the  RCS  vre.i ’'■'hti  ng  development). 

It  vra.s  also  possible  to  derive  the  ''ramer-Pao  bound  on  th.e  monopul 
variance.  This  bound  v/as  derived  and  is  ;iven  by: 

QT^  =  ^^V'2/N 

CM 

The  performance  of  the  :!cr>  v.'oi^iiting  technique  is  seen  to  uppi'oacii  t;ie  Cr.i.'..’a 
.Rao  bound  for  large  M. 

As  a  result,  tracking  filters  t.-ich  can  ifci;,.;it  Lte  injid.  ii.asuro- 
ments  should  use  a  weip.hting  which  is  proportional  -c  '..nc-  lecciVoa  signal 
pov.or, 

'i.2.3  Adaptive  Frequency  Selection,  irci;.  t:ie  ..c,„..c  lern.uiatio.i  •  a- 
'.'auar  glint,  as  discussed  above,  it  is  cxear  L.ial  t..o  a!!io.4nt  oi  gii  ls 
inversely  proportional  to  Uie  r  anar  cross  section  ‘..'f-o  o;  rne  scattcrer. 
'ihus ,  wlien  the  scutLenug  centers  on  i.nt  target  aie  n  pnase,  .yaeiUing  a 
maxi).'iuj!i  radar  return,  the  radar  balance  point  is  wiUiin  tl»e  physical  bounds 
oi  ^.he  .argbt,  whereas  when  tne  sc  •.tfenn  center  oui  f  phase,  yiei'im' 

r.i  r.iini  iuiT,  i-auar  return,  the  rauar  uaiauce  point  can  i'e  far  o  \  ..i  he 
&:.„jji.,il  hounas  C)'  tne  target.  .'.ne;'  the  U.U'get  ^.s  iru';/.-  rolai  i  -  to  t;.'- 
v;,  .veio  i;:th,  tiie  relative  pnasin/’,  oi  ti.e  ..cuLi.eriiii  centers  can  cnar..,c  evr 
wit.i  .^maii  cfiangOG  in  target  or-ientation,  so  tne  a:,  ount  oi  glint  can  nei  '  an 
upon  preciiv  nrget  orientation.  Sfnail  ctianges  in  rauar  iron  enc.v  e  ui 


affect  the  relative  phaein^j  of  the  returtia  from  the  various  soatteriiiiJ: 
centers,  thus  affectin^J  the  arnounc  of  glint. 


Because  of  the  frequency  dependence  of  the  relative  phasing  of 
scattering-center  returns,  it  would  appear  to  be  possible  to  maintain  a 
maximum  RCS,  and  thus  a  minimum  glint,  through  adaptive  modification  oi  i;nG 
radar  frequency.  This  possibility  has  been  recognized  for  a  long  time;  ..n 
fact  it  appears  in  one  of  the  first  [Dapers  devoted  to  the  analysis  of  radai' 
glint.*  A  review  of  miich  of  the  later  literature  on  radar  glinr,  indicator 
that  little  attention  seems  to  have  been  paid  this  ocssible  approach, 
although  a  great  deal  of  effort  has  been  devoted  to  t'^chniquos  in''^oivtn.'' 
random  frequency  variation  (or  nonrandom  variation  over  rang.es  large  r.r.ourh 
to  assure  decorrelation  of  the  radar  returns)  with  subsequent  averaging, 
thresholding,  filtering,  etc.  In  the  remainder  of  this  section  the  problem 
inherent  in  adaptive  frequency  selection,  a  problem  that  has  probably  .on- 
der  ’  attempts  at  developing  such  systems  futile,  is  di,-c'ass>^d. 

In  principle,  it  would  seem  re.,son£*ole  to  vary  frequency  on  a 
pulse-to-pulse  basis  until  a  relatively  large  HC'3  is  obtained.  (Hote  tnat 
a  truly  maximum  RCS  is  not  really  needed — what  is  important  is  to  avoid  nui 
in  RCS  that  lead  to  very  large  glint  values.)  Once  a  large  RCS  is  obta:. 
the  frequency  is  changed  in  small  steps,  of  the  order  of  i  to  4  Milt,  for 
typical  target  sizes,  with  the  direction  of  change  being  determined  by 
whether  the  RCS  is  increasing  or  decreasing.  Thus  if  a  clianrc  by  i  Miz 
(e.g.)  in  the  increasing-frequency  direction  loads  to  an  l  it  t'ease  lp.  b, 
then  the  frequency  is  increased  by  1  iHz  on  the  next  pulse,  ii  Pne  o  ae- 
creases,  the  frequency  is  reduced  by  2  MHz  on  the  next  pulse  to  cor!:>encate 
for  having  gone  the  wrong  way  witli  the  previous  pulse.  Ihis  form  of  adapPi 
frequency  selection,  a  simple  hili-olimbing  technique,  siiould  permn.  out:  to 
remain  necir  the  peak  RCS,  thus  reducing  the  amount  of  radar  glint.  The 
shifting  of  frequency  can  eventually  lead  to  the  Byslem  frequency  bounu, 


L.  Peters,  Jr.,  and  F.C.  Weimer,  "Tracking  radars  for  complex  targ,et 
Pr PC .  lEE ( London ) ,  vol.  110,  pp.  2149-21C2,  December  190,3. 


since  the  radar  cannot  operate  over  an  unrestricted  frequency  range.  Vihen 
such  a  bound  is  reached,  the  radar  must  again  use  the  random-frequency 
approach  to  acquire  a  new  peak  RCS. 

When  an  adaptive  technique  such  as  that  described  above  is  tried, 
it  is  found  that  the  pea!<  RCS  is  indeed  found  and  tracked  for  a  while  but 
then  reems  to  be  lost:  the  system  seems  to  track  small  RCS  values  rather  than 
large  ones.  Frequent  reacquisition  of  RCS  peaks  is  necessary.  The  reason 
for  this  problem  is  an  apparently  inherent  one  in  the  behavior  of  RCS  as  a 
function  of  aspect  angle  and  frequency.  If  one  thinks  of  RCS  as  an  elevation 
over  a  plane  in  which  coordinates  are  aspect  angle  and  radar  i^equency, 
then  what  one  seeks  is  the  frequency  that  produces  a  peak  RCS,  i.e.,  a  max¬ 
imum  elevation  on  the  surface  (and  hence  the  hill-climbing  algorithm  referred 
to  above).  If  the  maxima  of  this  surface  are  ridges,  then  as  aspect  angle 
changes  it  should  be  possible  to  vairy  the  frequency  in  a  manner  that  will 
keep  us  on  the  ridge  (assuming  the  PRI  is  short  enough  so  the  aspect-angle 
change  is  very  sirall).  If,  on  the  other  hand,  the  maxima  are  isolated  peaks, 
then  as  aspect  angle  changes  one  could  not  remain  on  a  peak  •  c.ianping 
frequency  in  the  nearly  continuous  manner  of  the  hill-climbing  technique 
described  above. 

To  permit  an  assessment  of  the  nature  of  the  RC'^  •'’irfare  -><=•  • 
function  of  aspect  angle  and  frequency,  a  simple  three-scn+rter  target  P' 
shown  in  Figure  11  was  chosen.  The  point  scatterers  have  scattering  ami^li- 
tudes  n.C,  and  D,  where  the  scattering  amplitude  is  the  na.uaro  root  of  RCb. 
The  discussion  is  limited  to  two  dimerjsions  since  the  addition  the  third 
dimension  adds  nothing  to  the  conclusions.  Initiull>  assicfo  scattering  am¬ 
plitude  B  =  0  so  there  aire  only  two  significant  scattering  centers.  The 
radair  cross  section  is  giver  by 

RCS  =  +  d2  +  2CPcos|^  .?k.it-d’coBo| 

where  k  is  the  wave  number  and  the  otlier  parameters  are  d  a'-',  ve 


Ficure  11.  THRFF-SCATTF RER  MODEL 
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To  remain  at  the  peak  of  the  RCS  function,  it  is  necessary  that 


2k(c+d)cos&  =  n2/T 

with  n  any  integer.  Expressing  k  in  the  form  2Tff/v,  where  f  is  the  radar 
fl'equency  and  v  is  the  velocity  of  propagation  (used  here  because  c  is  a 
length),  one  can  show  that  the  frequency  necessary  for  peak  RCS  is  given  by 

f  =  nv 

2(G+d)cosC 

There  are  n  possible  frequencies  because  of  the  multiple  peaks  of  the  cosine 
function.  For  c+d  =  5  meters,  the  values  of  f^^  for  selected  values  of  n  are 
shown  in  Figure  12.  For  a  given  value  of  n  the  variation  in  f  necessary  to 
stay  on  this  ridge  of  the  function  can  vary  widely  unless  the  aspect  angle  i 
restricted  to  lie  in  a  narrow  band.  If  a  frequency  limit  of  the  system  is 
reached,  it  is  necessary  to  shift  to  another  ridge  of  the  RCS  function,  i.e. 
to  a  different  value  of  n.  Still,  for  the  2-scatterer  case  it  seems  feasibl 
to  remain  on  such  a  ridge  with  only  periodic  shifts  to  new  values  of  n. 

It  is  obvious  from  the  result  given  above  that  the  intuitive 
hill-climbing  idea  for  frequency  adaptation  should  be  practical  for  the 
Lwo-scatcerer  target.  Practical  targets  often  consist  of  more  than  two  scat, 
tering  centers,  and  it  is  necessary  to  examine  the  behavior  of  the  RCS  conto 
for  such  targets.  Now  consider  the  3-scatterer  configuration  shown  in  1- i  ■  i..> 
11  with  B  made  nonzero.  The  RCS  is  given  by 

RCS  =  B'’  +  c2  p2  2RCcos  (2k  (  ccose  -  hsinO)j  t- 

2PDcos(21t(  dcosU  +  bslnO)j  tPCL'cos  (pk  ( c  vd  i  ■  tol 

now  one  has  three  terms  involving  cosine  functions,  and  it  seems  clear  tiiut 
one  can  no  longer  rentain  on  a  peak  of  Uie  Overall  function  oy  varying  the 
frequency,  since  there  are  three  teri.s  tnat  involve  the  asp  i,  onr.ir-  and 
tar.^et  dimensions  in  different  ways.  Thus  one  has  throe  sets  oa  ridge-anu- 
vailey  structures  that  are  superimposed,  and  the  locations  oi  peaks  depend  c 
the  I'jcations  on  the  frequency-angle  plane  at  idiich  three  riuges  superinij  e 
Figure  13  shows  the  ridge  orientations  for  the  ‘diree-componont  structure:  . 


f  -  GHz 

:'>CA'r‘TiP  Or.F  SHO^rNG  ridge  lines  and  80^+  peaks  (b=lm,  c=d'=2.5m) 


'ITie  set  of  curves  consisting  of  long  and  short  dashes  corresponds  to  scatter¬ 
ing  centers  C  and  D;  these  lines  are  exactly  the  same  ones  shown  on  the 

preceding  figiore  for  the  two  scatterer  model.  The  solid  curves  give  the 
ridge  locations  for  scatterers  B  and  C,  and  the  dashed  curves,  the  ridge  loca¬ 
tions  for  scatterers  3  and  D.  (.Vhere  the  ridges  coincide,  one  would  expect 
to  obtain  peaks  of  the  RCS  function.  For  the  set  of  scattering  amplitudes 
B-1,  «:.  =  2,  D=l,  shown  on  this  figure  are  approximate  contours  within 
ivhich  the  RCS  is  80%  or  more  of  its  maximur'  value. 

In  the  figure  above  all  of  the  ridges  lines  intersected  at  common 
points,  thus  producing  relatively  symmetrical  peaks.  This  intersection  of 
the  sets  of  ridges  is  a  result  of  assuming  c  =  d.  If  c  and  d  are  not  assurucd 

equal,  a  less  symmetric  set  of  peaks  occur  and  a  more  complicated  surface 

results.  Similar,  more-complicated  surfaces  v/ill  also  result  when  more  than 
three  scattering  centers  contribute. 

VJhen  the  presence  of  a  third  scattering  center  produces  isolated 
peaks  rataer  than  continuous  ridges,  one  can  see  that  the  peak -finding  (or 
hill-cliinbing)  technique  has  very  limited  applicability.  If  the  aspect  angle 
changes  by  enough  to  take  us  off  a  local  peal<,  the  adaptive  method  can  get 
trapped  in  a  low-RCS  region  for  many  pulses.  It  v/ill,  in  such  cases,  give 
significantly  poorer  perforinance  t?ian  a  random-frequency  technique  tliat  uaes 
thresholding  or  filtering.  The  only  way  the  adaptive  approach  car.  be  valid 
is  if  it  incorporates  a  means  of  recognizing  the  locations  of  peaks  anu 
sv'itching  to  a  random- frequency  tactic  when  a  peak  is  left,  thus  seeking  a 

1,'trv  ,  e  ak . 
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3.3  Simulation  Results 


3.3.1  Target  Scattering  Models.  To  permit  numerical  evaluation  of  radar 
glint  and  glint-reduction  techniques,  two  target  scattering  models  were 
developed.  The  first  model  was  for  an  F-15  fighter  aircraft.  Detailed 
dimensions  were  not  available,  so  figures  from  Janes'  were  used  in  conjunc¬ 
tion  with  measxa^ements  made  on  a  1/72  scale  model  (Revell  kit).  The  second 
scattering  model  was  for  a  cruise  missile.  Here  even  less  detailed  infor¬ 
mation  was  available  and  dimensions  were  based  upon  vdiat  information  could 
be  gleaned  from  Janes'  and  from  photographs  in  Aviation  VJeek  &  Space  Technol¬ 
ogy  articles;  the  TomahavAc  was  chosen  since  more  information  seemed  to  be 
available  for  it. 

For  both  cases,  assumptions  were  made  to  permit  reasonably  simple 
radar  scattering  analysis  to  be  used.  These  assumptions  included: 

•  Short  radar  wavelengths,  corresponding  to  frequencies  above 
3  GHz, 

•  Monostatic  radar  operation, 

•  Aspect  angles  limited  to  range  of  130°  in  elevation  and 
l50°  in  azimuth, 

•  Single  scattering  is  assumed,  i.e.,  no  multiple-diffraction 
effects  are  included, 

•  Shadowing  boundaries  were  crudely  approximated,  partially  tc 
simplify  the  analysis  and  partially  because  detailed  considera¬ 
tion  would  require  inclusion  of  multiple-diffraction  effects, 

•  Scattering  from  engine  ducts  was  modeled  in  a  very  crude  manner , 
partly  because  of  inadequate  knowledge  of  scattering  processes 
in  ducts  but  primarily  because  of  a  total  lack  of  information 
regarding  the  materials  and  geometry  of  the  engine  ducts 
involved  in  the  F-15  and  the  Tomahawk, 

More  detailed  descriptions  of  these  scattering  models  are  contained  in  the 
following  subsections. 


K-15  Scattering  Kodel.  An  outline  drawing  of  the  F-15  and  a 
list  of  its  major  dimensions  are  given  in  Figure  14  \Aich  has  been  taken  from 
Janes',  The  scattering  from  the  aircraft  was  modeled  by  81  scattering 
centers,  most  of  them  simple  point  scatterers  used  to  represent  scattering 
from  the  various  external  stores  on  the  aircraft.  Several  options  are 
available  in  the  computer  program  as  to  which  of  the  stores  are  included;  if 
ail  are  emitted,  the  aircraft  itself  is  mo'Jeled  by  31  scattering  centers.  A 
list  of  the  scattering  centers  is  given  in  Table  I;  some  of  the  repetitive 
items  are  omitted. 

The  manner  in  which  the  scattering  is  approximated  will  not  bo 
given  in  detail  here  because  of  the  complexity  and  length  of  discussion  tliat 
w'ould  be  needed.  Instead,  we  will  indicate  the  type  of  approximation  that 
is  used  for  each  oi'  the  classes  of  scattering  centers  on  the  aircraft. 

Edges  occur  at  engine-uucl  inlets,  leading  and  trailing  edges  of 
w-aj-,,  .»ud  tail  surfaces,  and  on  the  missile  support  structure.  Scattering 
.'rom  xhai'i.'  oiiges  is  approximated  through  use  of  tho  geometrical  theory  of 
diffraction,  which  permits  inclusion  of  polarization  effects,  for  a  tv;o- 
dimensional  edge  in  combination  with  a  pliysical-optics-type  approximation  to 
account  for  the  nonspecular  incidence  effects  on  a  finite-length  edge.  This 
type  of  scattering  approximation  should  be  reasonably  accurate  for  angles  not 
too  far  from  specular;  for  angles  that  tiro  far  from  specular,  the  errors  in¬ 
crease  but  the  amount  of  backseat  ter  decreases  so  that  tlie  errors  are  of 
•  jss  significance  tJian  would  otiierwiso  be  the  case. 

Scattering  fror:  the  scanning  antenna  in  tlie  nose  of  the  F-15  was 
approximated  by  tlie  scattering  from  a  flat  plate.  Scanning  of  the  antenna 
independent  of  aircraft  motion  is  possible.  Tliis  option  lias  generally  not 
been  us^d  in  tlie  computations  reported  fiere  (the  antenna  return  can  be 
switch  jd  ofi  for  any'  run). 

Scattering  from  the  junction  of  the  radomc  and  the  i'uselage  could, 
not  be  very  accurately  modeled  because  we  lacked  sufficient  information 


TABLE  1 


SCATTERING  CENTERS  ON  F-15  MODEL 
On  Aircraft 

Top  edge  of  port  engine  duct  intake 

Top  ed  .;e  of  starboard  engine  duct  intake 

Leading  edge,  port  v\/ing 

Leading  edge,  starboard  v;ing 

Trailing  edge  of  port  wing,  outer  section 

Trailing  edge  of  stai'board  wing,  outer  section 

Trailing  edge  of  port  wing,  inner  section 

Trailing  edge  of  starboard  v/ing,  inner  section 


Leading 

edge 

of 

horizontal 

stabilizer. 

inner 

section, 

port 

side 

Lead'!  ng 

edge 

of 

horizo’ital 

stabilizer, 

inner 

section, 

star bora d 

side 

be  :ding 

edge 

of 

horizontal 

stabilizer, 

outer 

section, 

port 

side 

.  ling 

edge 

of 

horizontal 

stabilizer, 

outer 

section, 

starboard 

side 

Trailing  edge  of  hori?:ontal  stagilizer,  port  side 
Trailing  edge  of  horizontal  s-'abilizer,  starboard  sid.e 
Leading  edge  of  port  vertical  stabilizer 
Lea.dlng  edge  of  starboard  vertical  stabilizer 
Trailing  edge  of  port  vertical  stabilizer 
Trailing  edge  of  starboard  vertical  stabilizer 
Scanning  antenna  in  nose  of  aircraft 
Junction  of  radome  and  fuselage,  port  side 
Ju'clion  of  radome  and  fuselage,  starboard  side 
d  j':i,;:lt 

Lower  edge  of  engine  duct  intaiie,  port  side 
Lov/er  edge  of  engine  duct  intake,  starboitrd  siue 
Scatterer  near  end  of  port  wing, 

^io  it'.crer  near  end  of  starboard  wing 
Scatterer  near  top  of  port  vertical  stabilizer 
Dcutlcrer  ncftr  top  of  starboard  vortioc.l  staV, i I i /.or 
Caiiniin  at  root  of  starboard  wing 
i  or-*  engine  di.K;t 
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TABLE  1 
Continued 


Starboard  engine  duct 


Stores 


Rear  edge  of  port  missile  mount 
Rear  edge  of  starboard  missile  mount 
Nose  of  Sparrow  missile  1 
First  set  of  fins,  Sparrow  missile  1 
Rear  set  of  fins,  Sparrow  missile  1 
Etc.  for  Sparrow  missiles  2,3,4 
Nose  of  Sidewinder  1 
First  set  of  fins,  Sidewinder  1 
Rear  set  of  fins,  Sidewinder  1 
Etc.  for  Sidewinders  2,3,  4 
Nose  of  bomb  1 
Tail  of  bomb  1 

Etc .  for  bombs  2  tlirough  12 
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about  the  aircraft,  '.v’e  therefore  used  point  scatterers  at  the  junction 
points  as  a  crude,  first-order  approximation  to  the  actioal  scattering. 

Scattering  from  the  cockpit  was  approximated  in  two  v;ays,  chosen 
by  an  input  option.  Either  the  use  of  a  conductive  canopy  could  be  assumed, 
in  which  case  there  is  no  scattering  from  either  the  interior  of  the  cockpit 
(because  it  is  shielded)  or  from  its  exterior  (because  we  are  restricted  to 
aspect  angles  for  which  specular  reflection  does  not  occur),  or  a  transpsu'ent 
canopy  could  be  assumed,  in  which  case  scattering  from  a  sphere  at  the  loca¬ 
tion  of  the  pilot's  head  is  assumed.  Additional  scattering  from  tiie  interior 
of  the  cockpit  also  probably  occurs  but  h.as  not  been  approximated  here. 

Additional  scattering  occurs  from  the  lower  edges  of  the  engine- 
duct  intakes;  either  sharp-edge  scattering  or  scattering  from  a  cylinder  is 
used  to  represent  this  scattering,  depending  on  the  elevation  angle  of  the 
111:,  of  sight.  Scattering  from  small  objects  at  the  ends  of  the  wings  and 
from  the  Ceinnon  on  the  starboard  wing  are  represented  by  point  scatterers. 

As  noted  above,  scattering  from  the  engine-ducts,  which  represents 
a  major  source  of  scattering  for  angles  near  nose-on,  could  not  be  very  ac¬ 
curately  modeled  because  of  a  lack  of  detailed  information  regarding  the 
ducts  as  well  as  fundamental  difficulties  in  modeling  duct  scattering  in 
any  case.  Basically,  we  assumed  a  fixed  scattering  amplitude  for  the  duct 
and  then  used  an  aspect-angle-dependent  factor  to  account  for  reduction  in 
illiimination  as  aspect  angle  varies.  It  is  possible  tliat  a  radar  absorbing 
grid  might  be  used  in  the  diet;  to  allow  for  this  case,  an  input  parameter 
can  be  set  to  reduce  the  duct  scattering  by  approximately  14  dB  to  approximate 
the  effects  of  reduced  duct  scattering. 

Scattering  from  the  missiles  and  bombs  that  can  be  carried  by  the 
aircraft  could  be  modeled  using  formulations  for  curved  surfaces  and  ogives 
in  combination  with  formulations  for  sharp  edges  such  as  are  used  for  edges 
on  the  aircraft.  Such  modeling  would  also  require  relatively  careful  inclu¬ 
sion  of  shadowing  effects  of  various  stores  by  each  otlier.  To  include  all 
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of  these  effects  would  require  a  very  complex  computer  program  which  would 
require  excessive  running  time.  Consequently,  we  have  adopted  a  much  simpler 
model:  simple  point  scatterers  are  used  for  each  set  of  fins  on  a  missile  or 
bomb  and  for  the  nose  of  each  missile  or  bomb.  While  of  limited  precision, 
this  form  of  approximation  should  be  reasonable  for  the  glint  analyses  of 
interest  here. 

The  key  characteristics  of  this  model  are  presented  in  Figure  15. 
Included  are  forward  aspect  angle  plots  of  the  target's  narrow  band  RCS  and 
the  first  and  second  central  moment  of  the  target's  RCS  distribution  in  each 
of  the  principle  measurement  dimensions.  These  moments  are  of  key  interest 
because  they  control  the  mean  and  variance  of  the  target  noise  in  the  cor¬ 
responding  radar  measurement  coordinates. 

The  data  presented  are  for  a  frequency  of  3  GHz,  It  is  of  interest 
to  note  that  the  aircraft's  RCS  is  fairly  insensitive  to  polarization.  This 
is  because  most  of  the  scattering  arises  from  engine  ducts  which  are  here 
modeled  as  polarization-insensitive  scatterers  rather  than  from  edges.  One 
would  generally  expect  polarization  diversity  to  be  ineffective  as  a  target 
noise  mitigation  technique. 

3. 3. 1.2  Cruise  Missile  Model.  Figure  16  is  a  sketch  of  the  simplified 
cruise  missile  configuration  that  was  modeled.  The  body  is  essentially  a 
cylinder  0.53m  in  diameter  and  5.49m  long.  The  nose  is  assumed  to  be  a 
(conducting  ellipsoid  of  revolution  with  major  axis  0.96m  and  minor  axes 
0.53m.  The  nature  of  the  duct  interior  was  not  know,  and  even  the  dimensions 
and  location  of  the  duct  entrance  had  to  be  inferred  ftom  various  photographs 
of  the  Tomahawk.  The  scattering  model  used  here  should  therefore  be  con¬ 
sidered  as  a  reasonable  representation  of  a  class  of  missle  configurations 
rather  than  as  an  accurate  model  for  the  Tomahawk . 

The  scattering  centers  used  to  represent  the  cruise  missile  are 
listed  in  Table  2.  Because  the  wings  and  tail  are  very  thin  and  have 
sharp  edges,  it  seems  reasonable  to  approximate  the  scattering  from  these 
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DBSM 


Length  =  .5.49m 
Wingspan  =  3.W7m 
Diameter  --  0.53ni 


GRirrSE  MISSIIE  USED  FOR  SCATTERING  MODEL 


edges  as  that  fyom  a  conducting  half  plane.  The  asymptotic  expression  for 
this  scattering,  using  the  geometrical  theory  of  diffraction,  is  the  sajne 
as  was  used  for  the  aircraft  aside  from  the  simplification  here  whereby  the 
wedge  angle  goes  to  zero  (i.e.,  the  half-plane  case).  For  the  leading 
edges  of  the  wings,  the  analysis  was  complicated  by  the  presence  of  rounded 
ends.  For  nonspecular  aspect  angles,  the  rounded  end  becomes  significant 
scattering  contributor;  although  the  scattering  amplitude  is  relatively 
small,  the  ends  of  -tdie  wings  are  the  scattering  centers  farthest  removed 
laterally  from  the  centroid  of  the  target  and  thus  can  produce  a  significant 
contribution  to  the  angular  glint  which  is  of  primary  significance  here. 
Analysis  of  scattering  from  a  straight  edge  with  a  curved  end  led  to  a  com¬ 
plex  expression  involving  a  uniform  asymptotic  expansion  of  the  scatter 5.ng 
integral;  the  uniform  asymptotic  expansion  is  necessary  here  because  the 
stationary-phase  point  of  the  integral  approaches  one  of  the  limits  of  th.e 
integration. 

The  nose  was  assumed  to  be  a  conducting  ellipsoid,  as  indicateii 
above.  Vie  did  not  have  information  on  the  location  of  the  radome,  but 
since  cruise  missiles  navigate  by  a  combination  of  inertial  navigation  and 
radio  altimeter  data,  it  seems  likely  that  the  antenna  is  directed  downwards 
and  thus  produces  little  scattering  in  tiie  forward-aspect  region.  Therefore, 
we  ignored  the  scattering  from  this  antenna. 

Scattering  from  the  interior  of  the  duct  v;as  treated  in  the  samt,- 
manner  as  the  duct  scattering  for  the  aircraft;  it  is  assumed  that  the 
scattering  is  reduced  proportional  to  the  intercepted  area  of  tVic  duct 
entrance.  Scattering  from  the  ring  discontinuity  at  the  duct  entrance  was 
assumed  equal  to  the  scattering  from  a  wire  ring  (i.e.,  an  eT.j..il  to  the 
intercepted  area  of  the  ring)  with  a  linear  dropoff  in  scattering  amplitude 
for  angles  off  nose-on;  this  rate  of  dropoff  depends  on  the  radar  wavelength, 
and  for  short  wavelengths  this  source  of  scattering  quickly  becomes  negligible 
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TABLE  2 


SCATTERING  CENTERS  ON  CRUISE-MISSILE  MODEL 

Leading  edge,  port  wing 
Leading  edge,  starboard  wing 
Trailing  edge,  port  wing 
Trailing  edge,  starboard  wing 
Leading  edge,  port  horizontal  tail 
Leading  edge,  starboard  horizontal  tail 
Trailing  edge,  port  horizontal  tail 
Trailing  edge,  starboard  horizontal  tail 
Leading  edge,  upper  ver-cical  tail 
Trailing  edge,  upper  vertical  tail 

Leading  edge,  lower  vertical  tail 

Trailing  edge,  lov/er  vertical  tail 

Nose 

Interior  of  duct 

Ring  discontinuity  at  duct  entrance 

3hadov/ing  of  the  scatterers  on  the  cruise-missile  niodel  is  very 
difficult  to  include  accurately,  llie  problem  is  that  most  of  the  soauterers 
are  edges  and  the  illuminated  portions  of  those  edges  are  extended  regions 
that  are  only  gradually  shadowed  as  the  aspect  angle  (e.g. )increases.  Because 
of  other  approximations  ti  at  had  to  be  made  in  the  model,  because  -  f  a  ia^i^ 
of  detailed  dimensional  information,  as  well  as  because  of  tlie  need  for  a 
model  simple  enougii  to  permit  very  rapid  coniputation  of  scattering  contr-  - 
butions,  we  adopted  a  very  simple  shadowing  argument.  Vie  assume  that  if  the 
aspect  angle  is  positive,  all  ed,;,e  sicatterers  on  the  starboard  side  of  the 
missile  are  shadowed,  and  vice  versa,  Tlius  at  0  degrees  aspect  angle,  i'-. 
abruptly  switch  from  the  scatterers  on  one  side  to  those  oti  the  otheu'.  A 
similar  rule  is  applied  for  scatterir.;;.  from  the  edges  on  the  upper  anu  thi^ 
lower  tail  surfaces:  for  negative  elevation  angles  the  lower-tail  edi^cs  are 
seen,  and  for  positive  aagies,  the  upper-tail  edges  are  seen,  'ibe  ’jn,  im 
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duct  is  assumed  to  be  seen  for  elevation  angles  less  tdian  +0.02r  and  to  be 
shadov/ed  for  elevation  angles  greater  than  +0.05r;  a  linear  dropoff  in 
scattering  for  angles  between  these  values  is  assumed.  Inclusion  of  the 
shadowing  effects  of  the  wings  on  the  duct  entrance  would  require  very  ex¬ 
tensive  computation  and  probably  would  not  significantly  affect  the  glint - 
analysis  results  in  any  case. 

The  key  characteristics  of  this  model  are  presented  in  Figures  17 
and  18.  Here  aspect  angle  plots  of  the  cruise  missile's  RCS,  RCS  centroid 
and  2*^^  central  RCS  moment  are  presented  for  both  horizontal  and  vertical 
polarizations.  Comparing  the  dual  polarization  RCS  plots  one  notes  a  imirked 
dissimilarity  in  behavior.  This  is  in  contrast  with  tlie  aircraft  results 
which  were  noted  to  be  largely  insensitive  to  polarization. 

While  thie  details  of  all  these  curves  will  not  be  discussed,  a 
brief  interpretation  of  the  behavior  of  the  cross-rani.e  RCS  centroid  will  be 
presented.  For  these  curves  the  missile  is  assumed  to  rotate  about  an  axis 
r\idway  between  nose  and  tail  and  the  centroid  location  is  measured  laterally 
boi.i  a  line  of  sight  passing  through  the  axis  of  rotation,  Tlius  if  all  of 
the  scattering  were  concentrated  at  the  axis  of  rotation  the  centroid  v;ould 
I’cmaiii  at  0 . 


Consider  first  the  horizontal-polarization  centroid.  There  is  a 
marked  trend  that  is  followed  by  the  centroid  for  most  aspect  angles,  Tliis 
trend  is  largely  a  result  of  the  gradual  motion  of  tho  duct  return,  although 
the  speculair  return  from  the  trailing  edge  of  the  vertical  tail  (which  has 

Scattering  for  horizontal  polarization)  also  contributes  to  t.he  glu  ft 
in  ctrntroid  location.  Actually,  the  return  from  the  tail  acts,  in  part,  .,o 
■  U'lzct  the  effect  (jf  scattering  from  tlio  nose,  wiiici:  would  in  turn  move  the 
centroid  towards  negative  values.  These  ttiree  principal  scatfcercrs  tb-uz 
produce  a  centroid  tliat  gradually  shifts  towards  positive-x  values.  Iioro 
■iCi-  a  large  ;)^-ak  ''xcursion  of  the  centroid  near  degT'ces.  This  shift 
.■  r  1*  ,:i'  s,  uoular  return  from  the  leading  edge  of  the  wing;  i|,i  j 

sc  .  tiiror  '  ’Comes  dominant  and  thus  shifts  the  centroid  towairis 
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tlie  leading  edge  of  the  wing.  The  smaller  petil^  near  11  degrees  is  similarly 
produced  by  the  leading  edge  of  the  horizontal  stabilizer.  As  this  edge  is 
sliorter,  its  RCS  is  smaller  and  the  centroid  excursion  is  smaller. 


For  vertical  polarization,  the  leading  edges  of  the  wing  ana  hori¬ 
zontal  fins  do  not  scatter,  so  the  pe-ilts  noted  in  tlie  horizontal-pol;ii'i,.ati 
case  do  not  appear.  Peal<ing  of  ti.e  return  at  0  degrees  Is  a  result  of  t;\o 
specular  returns  from  trailing  edges  of  v/lnp  and  tail  surfaces  v,tiich  are 
(assumed)  at  rigiic  angles  to  the  axis  of  the  missile.  Tlie  trend  o.oted  be¬ 
fore  is  no  iongei'  so  evii^letr. ,  because  the  trai litig-edge  specular  from  the 
tail  does  not  contribute  here;  instead,  wo  have  the  nose  ”ot:irn,  \idiic!.  lends 
to  move  th.e  nontr old.  towards  n'’gati\'o-x  values,  and  tlie  duet  rotia  :i,  v.tiieh 
ten-as  tc  move  it;  lowaras  the  positive-x  values. 


The  behavior  of  rn.a  radar  centroid  ttius  oaems  very  reasonable  in 
teri.is  of  the  suatlering  pr-ocesses  that  are  occurring  on  this  body.  Beca'.v 
of  the  approximations  inlierent  in  the  scattering  model  used  here,  it  would 
not  bo  reasonable  to  .at.to..nt  a  hi,;hly  detailed  analysis  of  centroid  motion. 
Our  aim  here  has  been  to  3U;.,^;est  that  the  behavior  tliat  has  been  observed  ws, 
iTi  fact,  consistent  wl'h  f'.-e  i.iodel  employed. 

c.J.l.d  (.tompai'isoii  wltli  HayieigJi  model,  oince  many  of  the  cand.'.ciite 
target  noise  reducticn  coace;;ts  iirc  basal  upon  the  resul  Ls  of  theoretical 
analysis;  which  assanieJ  a  Kayloijdi  ■•argot  model,  a  key  i  tom  ol'  interest  is 
to  /hat  extent  do  the  aircraft  and  oi'ulse  r.issilo  models  ap-proxiinate  a 
f.a/leigh  'tfirget.  To  effect  such  a  comp;irison,  histograms  of  the  RCf  val\  o  ■ 
oti'„ained  from  the  aircraft  and  cruise  missile  models  were  compared  with  a 
u,yleirli  distribution  and  tlie  results  presented  in  this  section. 

The  aircraft  R’Ch  ,  istograi.u;  in  Fi.-.urc  19  indicate  that  for  either 
I  olari:;a tion  tlic  f5c;ittering  arises  from  man.V  scattcrin>;  centers  and  thore- 
fci'o  r'ics  r.  Hayioig.!i  distrihation  'woll. 
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Unlilte  the  aircraft  model  vrtiich  showed  little  polarization  depen¬ 
dence,  the  cruise  missile  statistics  were  polarization  dependent.  The 
comparisons  with  a  Raylei^i  model  presented  in  Figure  20  illustrate  this. 

The  horizontal  polarization  data  appear  to  be  a  reasonable  fit  to  a  Raylei;^ 
model.  The  vertical  polarization  data,  however,  is  not  representative  of  a 
Rayleigh  distribution,  but  closely  resembles  a  two-scatterer  distribution. 
Since  the  vertical  polarization  RCS  is  strongly  dominated  by  the  cruise 
missile  nose  and  duct  returns  this  result  is  not  surprising. 

3.3.2  Simulation  Description.  In  section  3.3.1,  the  scattering  models 
that  were  developed  for  the  F-15  and  the  cruise  missile  were  described.  To 
permit  estimation  of  glint  errors  and  evaluation  of  glint-reduction  techniques 
it  was  necessary  to  develop  a  computer  program  that  uses  these  scattering 
models.  In  this  section  the  simulation  program  that  v/as  utilized  is  briefly 
described. 

Relatively  general  target  motions  were  designed  into  the  program. 

The  target  can  have  pitch,  roll,  and  yaw  angles  that  are  linearly  incremented 
with  time  (only  fixed-PRI  cases  are  treated)  and  can  have  arbitrary  starting 
points.  The  only  restriction  on  these  motions  is  that  the  line  of  sight 
must  produce  an  azimuth  angle  within  50°  of  nose-on  and  an  elevation  angle 
within  30°  of  nose-on.  This  limitation  resulted  from  tlie  fact  that  the 
targets  were  only  modeled  over  this  region. 

For  any  run,  three  frequencies  are  used  for  each  aspect  angle;  in 
this  way,  it  is  possible  to  evaluate  the  effects  of  various  frequency  be¬ 
haviors  on  glint  with  all  other  parameters  held  fixed.  The  first  set  of 
outputs  is  for  a  fixed  frequency,  i.e.,  one  which  is  exactly  the  same  for 
every  pulse.  The  second  set  of  outputs  is  for  a  randomly  varying  fre¬ 
quency.  Here  the  frequency  is  changed  from  pulse  to  pulse  on  a  pseudo¬ 
random  basis:  there  is  a  uniform  probability  the  frequency  has  any  value 
between  a  specified  minimum  and  a  specified  maximum  frequency.  The  minimum 
and  maximiun  frequencies  need  not  include  the  fixed  frequency  discussed 


666 


Lgure  2( 


above.  To  make  the  random-frequency  case  more  realistic,  a  frequency- 
quantization  value  is  also  assigned.  For  example,  a  frequency  range  from 
9  to  10  GHz  in  1-MHz  steps  may  be  specified.  There  would  then  be  1001  pos¬ 
sible  frequencies  and  they  would  be  selected  on  a  uniformly  distributed 
random  variable  basis. 

The  third  frequency  is  determined  using  an  adaptive  selection 
technique.  Here  the  maximum  and  minimum  frequencies  and  the  frequency  step 
are  the  same  as  were  assigned  for  the  random-frequencies  case,  but  an  attempt 
is  made  to  optimize  the  frequency  on  the  basis  of  the  observed  RCS  of  the 
target.  Recall  (Section  3.2.3)  that  adaptive  frequency  selection  is  expected 
to  be  difficult  because  of  the  isolated  nature  of  the  RCS  peaks  as  a  function 
of  frequency  and  aspect  angle.  Decause  of  limited  resources  and  some  skep¬ 
ticism  concerning  the  ultirfiate  utility  of  ft'equency  adaptation,  only  two 
adaptive  techniques  were  tested.  Of  these  two  techniques,  which  are  described 
below,  the  second  technique  performed  better. 

The  first  of  the  adaptive  frequency-selection  techniques  that  was 
tried  uses  random  frequency  selection  for  some  specified  (e.g,,  5)  number  of 
pulses  and  then  chooses  the  frequency  corresponding  to  the  largest  RCS.  nie 
frequency  is  then  stepped  (by  the  specified  step  size)  and  the  RCS  compared. 

If  the  RCS  has  increased,  the  frequency  is  stepped  by  the  same  amount  in  the 
same  direction.  If  tlie  RCS  has  decreased,  the  frequency  is  stepped  in  the 
opposite  direction  by  tv/ice  the  frequency-step  size.  Tiiis  procedure  continues, 
attempting  to  remain  near  tlie  peak  of  the  RCS  function.  As  shown  earlier, 
when  there  are  more  than  two  scattering  centers  on  the  target,  the  peaks 
can  become  isolated  rather  than  being  ridges  along  which  a  hill-climbing 
algorithm  can  be  expected  to  maintain  a  high  RCS  value.  Because  the  target 
as/icct  angle  variations  are  generally  not  controllable,  one  can  be  forced  off 
a  peal-:  in  spite  of  the  small-step  frequency  variations.  To  avoid  this  problem, 
whereby  one  could  be  trapped  in  a  lov/-RCS  region,  every  N'th  pulse  a  randomly 
chosen  frequency  (N^is  an  input  p<u'ameter  of  the  order  of  5  or  10)  is  used. 

If  the  RCS  at  the  randomly  chosen  frequency  is  greater  than  the  latest  RCS 
obtained  using  the  adaptive  frequency  selection  value,  the  adaptive  frequency 
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is  sv;itched  so  that  it  begins  to  v/ork  around  a  nevi  value.  In  this  way,  if 
one  is  getting  forced  into  a  "valley"  region  there  is  a  chance  that  a  randoi.i 
frequency  will  produce  an  RCS  nearer  a  peak.  The  nain  difficulty  v;ith  this 
method  is  that  in  some  cases  the  adaptive  method  gets  into  a  "valley"  region 
and  the  randomly  cliosen  frequency  also  is  in  a  lov/-RCS  region.  Many  pulses 
may  go  by  before  one  gets  back  to  a  hlgh-RCS  region.  Another  problem  is  that 
the  use  of  the  random  frequency,  when  vie  are  able  to  trad:  the  desired  fre¬ 
quency  adaptively,  leads  to  a  loss  of  usable  data. 

The  second  type  of  adaptive  processing  uses  the  same  principle 
of  stepping  the  frequency  toi;ards  larger  RCS  values.  Instead  of  trying  a 
random  frequency  every  H  pulses,  the  adaptive  method  is  used  until  the  RCS 
falls  to  a  level  1/3  that  of  the  largest  RCS  to  date.  Actually,  this 
largest  RCS  -value  is  allowed  to  decay,  to  prevent  trapping  of  -the  system  by 
an  abnormally  large  RCS  value  (e.g. ,  a  specular  return  only  seen  once),  using 
a  digital  filter.  If  the  observed  RCS  value  is  less  than  1/3  the  peak  (ac¬ 
tually,  filtered-peak )  value,  random  frequencies  are  used  until  the  RCS  is 
large  enough  for  the  adaptive  processing  to  be  resumed.  In  -this  way  one 
avoids  losing  good  data  -through  a  requirement  to  use  random  frequencies 
every  N'th  pulse,  but  now  requires  the  use  of  the  random-ft'equency  me-thod  for 
as  many  pulses  as  necessary  to  get  back  to  a  peak.  Hiis  second  type  of 
adaptive  processing  worked  quite  well  and  is  used  in  the  numerical  examples 
presented  later. 

For  each  pulse,  the  orientation  of  -the  -target  is  computed,  along 
with  -the  coordinate  -transformation  matrices  that  are  needed.  The  -three 
frequencies  to  be  used  are  then  computed  using  -the  techniques  described 
above.  For  each  fi?equency,  -the  RCS  centroid  and  -the  second  moment  of  the 
radar  centroid  is  computed.  Radar  glint  is  found  at  -the  same  time.  The 
glint  values  can  be  found  using  ei-ther  of  two  options:  it  can  be  found 
relative  to  -the  fixed  geometrical  center  of  -the  -target  (a  point  on  -the  cen¬ 
terline  of  the  target  and  a  specified  distance,  usually  '4  -the  length,  behind 
the  nose),  or  it  can  be  found  relative  to  the  target  centroid.  The  cen-troid 
and  glint  values  are  always  found  in  three  directions:  x,z,  and  range.  These 
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distances  are  measured  relative  to  the  radar  coordinates  with  x  a  horizontal 
displacement,  z  a  vertical  displacement,  Eind  range  measured  along  the  line 
of  si^t,  vrtiich  is  assumed  to  pass  through  the  previously  defined  geometrical 
center  of  the  target.  In  addition  to  the  centroid  and  glint  values,  sums  of 
these  quantities  and  their  squares  (in  some  cases)  are  stored  for  later  use 
in  computing  means  and  standard  deviations. 

An  option  permits  printing  out  pulse-by-pulse  data  in  one  of  two 
formats.  In  either  case,  the  pulse  number,  azimuth  angle,  and  elevation 
angle  are  printed.  For  each  frequency  the  RCS  is  also  printed.  For  the  first 
option,  the  x,  z,  and  range  glint  values  are  printed  for  each  frequency; 
either  the  glint  relative  to  the  geometrical  center  of  the  target  or  the  glint 
relative  to  the  radar  centroid  can  be  printed.  After  the  last  pulse  has  been 
processed,  the  table  concludes  with  the  mean  RCS  level  for  each  of  the  three 
frequency  sets.  This  mean  RCS  is  useful  for  evaluting  the  effectiveness  of 
adaptive-fhequency  processing,  as  a  highly  successful  adaptive  choice  of 
frequency  would  lead  to  an  increased  mean-RCS  level. 

The  second  option  for  the  pulse-by-pulse  table  again  prints  aspect 
angle,  frequency,  and  RCS  information  just  as  for  the  first  option.  In  the 
x,z,  and  range  columns,  instead  of  glint  values  the  locations  of  the  cen¬ 
troids  are  printed.  In  a  second  line,  below  the  centroids,  the  quantity 
(second  moment  -  centroid2)0*5  is  printed;  this  quantity  is  a  measure  of  the 
lateral  (or,  in  the  case  of  range,  longitudinal)  spread  of  the  scattering 
centers  on  the  target.  The  glint  computations  are  unchanged  by  this  choice 
of  pulse-by-pulse  output,  and  the  subsequent  summary  table  will  be  the  same 
for  this  case  as  for  the  first  choice  of  pulse-by-pulse  output. 

Statistical  information  on  the  target  scattering  characteristics 
is  given  in  a  table  that  follows  the  pulse-by-pulse  table  (when  it  is  used). 
First,  the  program  gives  data  relating  to  the  radar  centroid  of  the  target. 
These  data  are  given  for  each  of  the  three  frequency  sets  (i.e.,  fixed,  ran¬ 
domly  varied,  and  adaptively  varied)  and  for  each  of  the  three  corrdinates 
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(x,  z,  and  range).  Two  quantities  are  printed  for  each  centroid:  the  mean, 
and  the  second  central  moment. 

The  remainder  of  the  summary  table  contains  statistics  on  the  glint 
assuming  various  methods  of  processing  the  glint  information.  All  of  the 
glint  values  used  in  the  statistics  are  computed  relative  to  either  of 
two  references,  as  mentioned  above,  the  geometrical  center  of  the  target, 
which  remains  fixed  in  location,  or  the  radar  centroid,  vhich  slowly  varies 
in  position  on  the  target  as  aspect  angle  varies.  A  fixed  threshold  level, 
e.g, ,  -20  dBsm,  is  used  for  all  of  the  glint  processing.  The  reason  for  this 
threshold  is  that  in  practice  a  radar  has  limited  sensitivity  and  thus  does 
not  obtain  data  when  the  target  RCS  is  below  some  minimum  level;  for  very  low 
RCS  values,  the  glint  takes  on  very  large  values  that  would  not  be  seen  in 
practice.  This  fixed  threshold  level  can  be  made  very  low  if  inclusion  of 
all  glint  values  should  be  desired. 

Using  only  the  fixed  RCS  threshold,  for  each  frequency  and  for  each 
1  -hree  coordinates,  the  minimum  and  maximum  glint  values,  the  mean 
:  lint,  and  the  standard  deviation  of  the  glint  values  are  printed.  These 
results  Indicate  the  amount  of  glint  under  the  assumption  that  no  special 
processing  to  reduce  glint  is  performed. 

Next  the  glint  statistics  for  three  types  of  weighted  averages 

are  computed.  In  each  case  a  running  n-point  average  is  formed  and  the  mean 

and  standard  deviation  of  the  resulting  averaged  values  are  calculated.  This 

is  uone  in  each  of  the  three  measurement  coordinates  for  each  of  the  three 

frequency  selection  schemes.  The  three  weighting  techniques  implemented  are: 

uniform  wei^^iting;  RCS  wei^ting  which  was  shown  in  Section  3.2.2  to  be 

optirniim  for  Rayleigh  targets;  and  a  technique  termed  SLA  (Selection  of  the 

largest  Amplitude)  in  vhich  only  the  glint  measurement  corresponding  to  the 

largest  received  signal  is  selected.  The  SLA  technique  has  been  reported  to 

be  optimum  under  some  conditions.*  For  each  computer  run  up  to  5  different 

values  of  n  may  he  specified  by  the  user. _ 

*  J.K.  Loomis  and  E.R.  Graf,  "Frequency-Agile  Processing  to  Reduce  Radar 
Uint  pointing  Error,"  Trans.  IEEE,  vol.  AES-10,  Nov.  1974,  pp.  811-820 
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For  sitviatlons  in  vftilch  a  radar  tracking  filter  is  not  designed  to 
individually  wei^t  the  raw  radar  measurements,  data  editing  of  the  basis  of 
RCS  thresholding  may  be  used  to  reduce  glint  errors.  To  assess  the  effec¬ 
tiveness  of  such  techniques,  as  well  as  verifying  the  theoretical  results  of 
Section  3.2.2.,  the  simulation  program  also  computes  glint  statistics  for  5 
RCS  threshold  settings.  As  in  the  other  cases,  both  the  resulting  mean  and 
standard  deviation  are  computed  for  each  of  the  three  measurement  corrdinates 
and  each  of  the  three  frequency  selection  schemes.  In  addition  statistics  on 
the  percent  of  the  data  accepted  at  each  threshold  value  is  recorded.  This 
data  is  required  to  normalize  the  resulting  statistics  by  the  probability  of 
accepting  measurements  at  a  given  threshold  level. 

3.3.3  Simulation  Statistics.  This  section  contains  key  simulation 
results  for  the  two  specific  target  models  developed  and  described  in 
Section  3.3.1.  Generally  speaking,  when  several  dominant  scattering  centers 
existed,  the  target  scattering  was  approximately  Rayleigh  and  all  the  theo¬ 
retical  results  for  a  Rayleigh  tairget  developed  in  Sections  3.2.1  and  3.2.2 
were  found  to  apply.  In  these  situations  random  frequency  selection  and  RCS 
weighted  averages  produced  the  greatest  reduction  in  target  noise. 

3. 3. 3.1  Aircraft  Statistics.  The  simulation  was  used  to  gather  aircraft 
target  noise  statistics  at  both  horizontal  and  vertical  polarizations  at  a 
radar  frequency  of  approximately  3  (5Iz.  The  results  for  both  polarizations 
were  nearly  identical  so  only  specific  results  for  horizontal  polarization 
will  be  discussed. 

For  all  runs  a  -20  dBsm  "detection  threshold"  was  employed.  That 
is,  only  data  associated  with  aircraft  cross  sections  in  excess  of  this  value 
were  included  in  the  statistics.  For  the  data  presented,  the  mean  aiixiraft 
RCS  was  2.25  dBsm  i^ich  resulted  in  a  normalized  threshold  value  (i.e., 
threshold  value+mean  RCS)  of  x  =  .006.  According  to  the  theoretical  results 
of  Section  3.2.2  the  normalized  single  pulse  target  noise  varlauice  should  be 


This  value  is  noted  for  reference  on  Figure  21  which  presents  the  normalized 
target  noise  in  the  range  dimension  as  a  function  of  the  number  of  pulses 
averaged  for  each  of  the  three  frequency  selection  schemes.  In  all  cases 
RCS  weighting  was  employed.  From  the  curves  it  is  evident  that  of  the  three 
frequency  selection  schemes,  the  fixed  frequency  option  was  noticeably  the 
poorest.  Both  the  random  and  adaptive  frequency  options  produced  results 
near  the  Cramer-Rao  bound  of  1/N,  with  tiie  random  technique  performing 
slightly  better. 

For  random  frequency  selection.  Figure  22  presents  a  comparison 
of  averaging  techniques  for  the  same  simulation  run.  Here  uniform  weighting 
performed  the  worst  and  RCS  wei^ting  the  best.  Tliis  result  is  not  sur¬ 
prising  since  RCS  weighting  is  optimum  for  a  Raylei^  target  and  the  aircraft 
statistics  were  shown  to  be  nearly  Rayleigh  in  Section  3. 3. 1.3.  Similar 
results  were  also  obtained  in  the  other  measurement  dimensions  and  also  for 
the  other  frequency  selection  schemes. 


Figure  21  EFFECTS  OF  FREQUENCY  SELECTION  ON  F-15  MODEL  DATA 


NUMBER  OF  PULSES  AVERAGED 

Figure  22  EFFECTS  OF  WEIGHTING  ON  F-15  MODEL  DATA 

Finally,  results  dealing  with  the  impact  of  data  editing  on  the 

smoothed  position  of  a  target  under  track  are  presented  in  Figure  23,  Ihe 

quantities  graphed  here  are  the  normalized  target  noise  variance  in 
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each  of  the  three  measurement  dimensions,  further  normalized  by  the  prob¬ 
ability  of  obtaining  a  measurement  above  the  indicated  threshold  value.  The 
rationale  for  this  statistic  comes  from  the  fact  that  the  smoothed  target 
position  variance  varies  directly  with  the  target  noise  variance  and 
inversely  with  the  expected  number  of  measurements  used  in  the  filtering 
process,  NPp,  Generally  as  the  threshold  value  x  is  increased,  the  target 
noise  variance  is  reduced  but  so  is  the  number  of  measurements  exceeding 
the  threshold.  The  theoretical  results  of  Section  3.2.2  indicated  that  for 
a  Rayleigh  target  (see  Figure  10)  the  statistic  of  interest  has  a  minimum 
in  the  region  of  x  =  .6.  The  simulation  results  are  observed  to  generally 
support  this  prediction.  For  a  Raylei^  target,  the  probability  of  detection 
given  a  threshold  value  of  x  is  EXP(-x),  hence  for  x  =  .6  only  about  55% 
of  the  raw  radeir  measurements  should  be  used  for  tracking  with  algorithms 
which  do  not  permit  input  data  wei^jiting. 
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Horizontal  Polarization 
Random  Frequency  Selection 


NORMALIZED  THRESHOLD  -  x 

Figure  23  EFFECTS  OF  THRESHOLDING  ON  F-15  MODEL  DATA 


In  addition  to  the  frequency  selection  techniques  cited  above,  a 
limited  simulation  analysis  concerning  the  utility  of  polarization  diversity 
was  also  conducted.  Specifically  the  simulation  was  rian  for  both  horizontal 
and  vertical  polarization  at  the  same  frequency  and  the  pulse-by-pulse  meas¬ 
urements  stored.  In  addition,  a  third  run  at  horizontal  polarization  was 
made  at  a  frequency  offset  by  5  MHz  and  the  pulse-by-pulse  data  also  stored. 


m 


With  this  set  of  data,  RCS  weighted-tvro  pulse  averages  were  computed  for  both 
dual  polarization  and  dual  fl'equency.  Ihe  results,  which  are  presented  in 
Table  3,  clearly  indicate  that  dual  frequency  diversity  was  superior.  This 
was  not  an  unexpected  result  for  the  aircraft  model  because  of  the  hi^  degree 
of  correlation  (i.e.  over  80%)  between  the  dual  polarization  RCS.  On  the 
other  hand  the  5  MHz  change  in  frequency  was  sufficient  to  largely  decorrelate 
the  target  RCS  and  therefore  increase  the  performance  of  the  dual  diversity 
scheme . 


Table  3  Diversity  Comparisons 


^ Dual  Frequency 
0" Dual  Polarization 


Range 


.66 


Azimuth 


.70 


Elevation 


.67 


3. 3. 3. 2  Cruise  Missile  Statistics.  Using  the  cruise  missile  model  in 
the  simulation, analysis  similar  to  that  conducted  on  the  aircraft  model  was 
performed,  but  only  at  a  nominal  frequency  of  9.5  GHz,  These  results  are 
highlighted  in  this  section. 

The  first  general  result  pertains  to  the  pulse-to-pulse  frequency 
selection  procedure.  As  in  the  aircraft  runs,  compared  were  fixed  frequency, 
adaptive  frequency  and  random  frequency  selection  techniques.  Without 
expection  the  random  frequency  selection  consistently  provided  the  best 
performance  vrfiile  the  fixed  frequency  performed  the  worst.  Figure  24  presents 
a  representative  example  of  this  result.  Here,  for  horizontal  polarization, 
the  normalized  target  noise  variances  in  the  range  dimension  are  compared  for 
the  three  frequency  selection  techniques.  The  results  for  the  random  fre¬ 
quency  selection  are  not  far  from  the  Cramer-Rao  bound  for  Rayiei^  targets. 

The  next  result  of  interest  pertains  to  the  preferred  weighting 
technique  for  averaging  measurements.  In  the  theoretical  results  of  Section 
3.2.2  it  was  shown  that  RCS  weighting  was  optimum  for  a  Rayleigh  target. 
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Figure  24  EFFECTS  OF  FREQUENCY  SELECTION  ON  CRUISE  MISSILE  DATA 


Section  3. 3. 1.3  subsequently  showed  the  cruise  missile  statistics  to  be  es¬ 
sentially  Rayleif^  for  horizontal  polarization  but  not  for  vertical  polariza¬ 
tion  in  the  0°  elevation  cut  analyzed.  Representative  results  for  the  cruise 
missile  are  shown  in  Figure  25  and  26.  For  horizontal  polarization,  RCS 
weighting  is  shown  to  be  best,  which  is  consistent  with  the  Rayleigh  model. 
For  vertical  polarization,  however,  RCS  weighting  did  not  produce  the  best 
result.  Recall  that  for  this  polarization  the  target  RCS  statistics  were 
dominated  by  two  scattering  centers,  namely  the  missile  nose  and  engine 
intake.  For  this  case  a  stronger  RCS  weighting  is  evidently  optimum.  The 
selection  of  the  measurement  associated  with  the  largest  received  sunplitude 
(SLA)  was  previously  shoim  to  perform  best  on  simple  targets,*  and  indeed  it 


*  J.M,  Loomis  and  E.R,  Graf,  "Frequency-Agile  Processing  to  Reduce  Radar 
Glint  pointing  Error",  Trans.  IEEE,  vol.  AES-10,  Nov,  1974,  pp.  811-820 
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Figure  25  EFFECTS  OF  WEIGHTING  ON  HH  CRUISE  MISSILE  DATA 


was  found  to  perform  better  than  RCS  weighting  for  the  case  at  hand. 

Representative  results  of  the  effects  of  threshold  setting  on  data 
editing  for  tracking  are  presented  in  Figure  27.  In  general  the  best  thresh¬ 
old  setting  for  horizontal  polarization  was  near  the  theoretical  threshold 
value  of  .6  and  produced  a  normalized  variance  very  near  the  predicted  value 
of  1.5.  Since  the  vertical  polarization  case  was  far  from  a  Rayleigh  situ¬ 
ation,  it  was  not  surprising  that  a  similar  optimum  threshold  did  not  exist. 
In  keeping  with  the  SLA  finding  of  Figure  26,  the  threshold  is  probably  the 
one  that  yields  the  single  SLA  value. 
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Figure  26  EFFECTS  OF  WEIGHTING  ON  W  CRUISE  MISSILE  DATA 


The  final  result  of  interest  pertains  to  a  comparison  of  dual 
polarization  and  frequency  diversity  as  a  target  noise  reduction  technique. 
Unlike  the  aircraft  situation,  the  dual-polarization  cruise-missile  cross 
sections  were  not  highly  coxnrelated  and  polarization  diversity  proved  nearly 
as  effective  as  frequency  diversity. 
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3.4  Operational  Implications.  In  general,  target  noise  effects  mus . 
considered  in  the  radar  measurement  error  budget  i4ien  this  error  component 
becomes  comparable  to  the  thermal  noise  component.  In  Section  3.2.2,  w.- 
combined  glint  error  variance  for  a  Rayleigh  target  was  shou'  to  be  x 

imatelv 


a 


2 


.103  ©B,/ 

SNR  ■*'  ^ 


Here  the  first  term  represents  the  thermal  noise  component,  and  the  second 
term,  the  target  nose  component.  Equating  these  two  terms  and  graphing  the 
resulting  equality  one  obtains  Figure  28.  Here  the  ratio  of  the  target's 
angular  extent  parameter  Wq,  normalized  by  the  radar  beamwidth  is  plotted 
vs  signal-to-noise  ratio.  For  situations  above  the  curve,  target  glint 
effects  must  be  considered.  Completely  analogous  results  also  exist  for 
target  range  noise. 


Figure  28  CONTOUR  OF  EQUAL  THERMAL  &  TARGET  NOISE 


In  terms  of  diversity  schemes  for  decon^elating  pulse-to-puise 
target  noise,  frequency  diversity  appears  to  be  the  best  choice.  Generally 
speaking,  temporal  diversity  doesn't  produce  sufficient  decorrelation  and 
depending  upon  the  scattering  characteristics  of  a  target,  polarization 
diversity  may  not  be  effective  either.  As  far  as  fi'equency  selection  schemes 
go,  pseudo-random  frequency  selection  appeared  to  be  the  best  all  aroiand 
technique. 

Given  the  means  of  obtaining  uncorrelated  pulse-to-pulse  target 

noise  retvirns,  the  question  of  data  editing  or  averaging  must  be  considered. 

For  complex  targets  such  as  aircraft,  the  target's  RCS  statistics  can  be 

assumed  to  be  nearly  Reiyleigii,  and  the  theoretical  Rayleigh  results  utilized. 

For  tracking  filters  which  do  not  permit  weighting  of  the  input  measurement 

data,  a  data  editing  scheme  on  the  basis  of  the  received  target  power  was 

shown  very  effective.  Here  a  threshold  setting  equal  to  six  tenths  of  the 

average  target  power  is  recommended.  The  average  variance  of  the  radar 

measurements  accepted  in  this  manner  was  shown  to  be  approximately  .82  w|/2, 

where  is  the  target  extent  parameter  in  the  radar  measurement  dimension  of 

interest.  Since  only  about  55%  of  the  available  measurements  are  accepted  by 

the  threshold  test,  the  resultant  smoothed  target  position  variance  is 

V/2 

proportional  to  1.5  vrtiere  N  is  the  total  number  of  unthresholded  measure¬ 
ments  . 

For  those  situations  in  vniich  input  data  weighting  is  permissible, 
weighting  on  the  basis  of  the  signal  power  was  shown  optimum  for  Rayleigh 
targets.  This  type  of  wei^ting  produced  "averaged"  target  variamces  near 
the  Cramer-Rao  bound  of 

For  targets  vhose  scattering  characteristics  were  more  representative 
of  a  deterministic  two  scatterer  target  than  a  Raylei^  target,  heavier  RCS 
data  editing  and  weighting  is  apparently  optimum.  In  this  situation  better 
performance  was  obtained  with  the  SLA  technique  than  with  RCS  weighting. 
Effectiveness  of  intermediate  techniques  were  not  evaluated.  As  a  result, 
while  RCS  weighting  can  still  be  utilized  quite  effectively,  better  techniques 
do  exist  for  non-Rayleigh  targets. 


4.  KEY  FINDINGS  AND  RECOMMENDATIONS 


Addressed  in  this  study  were  aspects  of  target  noise  effects  as 
tliey  pertain  to  the  long  range  detection  problem  and  the  short  range  very 
precise  track  problem.  By  far,  the  major  effort  and  most  significant 
findings  dealt  with  the  precision  track  problem. 

For  the  analytic  studies  in  support  of  both  of  these  problems,  the 
target  amplitude  statistics  were  assumed  Rayleigh.  Generally  4  to  6  near 
equal  target  scatterers  were  found  adequate  for  near  Rayleigh  statistics.  A 
fairly  detailed  modeling  of  an  F-15  aircraft  also  produced  an  essentially 
Raylei^  target  behavior.  The  detection  analysis  j)erformed  was  limited  to  a 
search  radar  that  noncoherently  Integrated  target  returns.  Generally  speak¬ 
ing,  detection  performance  was  shown  to  degrade  with  increasing  pulse-to- 
pulse  correlation.  While  the  general  comparison  between  unity  pulse-to-puise 
correlation  (Swerling  I)  and  zero  pulse-to-pulse  correlation  (Swerling  Jt) 
is  well  documented  the  results  for  intermediate  correlation  values  ar'^ 
b.'i  Laved  new  but  not  profound. 

Several  new  and  novel  results  were  obtained,  hov/ever,  for  the 
targot  noise  problem  associated  with  very  precise  radar  measurements.  Hiesf. 
r. 'Suits  pertain  to  the  radar  measurement  errors  of  glint  and  the  analog, ous 
range  noise  error.  In  both  instances,  the  key  target  parameters  of  intent  be 
are  the  target  RCS  centroids^,  in  each  radar  measurement  coordinate,  ujid 
the  second  central  RCS  moments  W^,  in  these  same  coordinates.  For  high  SKu 
co.idit  ions,  the  average  target  measurement  location  was  shown  to  be  U;t; 
corresponding  location  of  the  target  RCS  centroid.  Tho  Cramer-Rao  bound  on 
tiie  associated  target  noise  variance  was  shown  to  be  W^/2N,  where  N  is  th.' 
nuiaber  of  independent  measurements  averaged.  When  averaging  independent 
measurements  under  high  SNR  conditions,  RCS  weighting  was  shown  to  be  opti-.u 
a. id  resulted  in  a  variance  of  w2/2(N-l). 


Finally,  in  those  instances  in  which  wei^tlng  is  impractical, 
data  editing  on  the  basis  of  an  RCS  threshold  was  shown  to  be  very  effective 
in  reducing  target  noise  variances.  A  threshold  of  six  tenths  of  the  mean 
target  RCS  was  shown  to  be  optimum. 

In  general  all  these  Rayleigh  results  were  supported  via  simulation 
analyses  using  the  aircraft  model.  Similar  results  for  a  cruise  missile 
model  indicate  that  a  heavier  RCS  wei^ting  is  optimum  for  simpler  non- 
Rayleigh  targets,  i.e,  two  point  scatterers.  This  result  was  not  pursued  to 
any  great  extent  and  may  warrant  additional  investigations  in  the  furture. 

To  provide  uncorrelated  pulse-to-pulse  target  noise  statistics, 
frequency  diversity  was  found  to  be  preferable  to  both  temporal  and 
polarization  diversity.  Pseudo  random  frequency  selection  was  Shown  to 
be  more  effective  in  reducing  target  noise  errors  than  either  fixed  frequency 
techniques  or  a  class  of  adaptive  frequency  selection  techniques  investigat’o. 
An  analysis  of  the  basic  problem  of  trying  to  adapt  the  radar  frequency  in 
sich  a  way  as  to  maximize  the  target  RCS  and  hence  minimize  target  glint, 
raised  serious  questions  as  to  the  feasibility  of  such  a  concept. 


APPENDIX  A.l  DERIVATION  OF  CORRELATION  COEFFICIENT  FOR  ROTATING  BODY, 


In  Section  2.1.2  an  expression  for  the  correlation  between  two  returns 
at  different  frequencies  was  developed.  A  similar  expression  can  be  devel¬ 
oped  for  two  returns  from  a  rotating  body  taken  at  different  times.  The 
geometry  is  as  shown  in  Figure  2  with  the  body  rotating  about  the  origin  of 
the  Cartesian  system. 

With  the  frequency  fixed,  the  complex  narrowband  return (assuming  small 
rotations)  is  represented  as  a  function  of  time: 


A(t)  =  E:<p(-j47rR 


/2W 

r  ^  /  y 


0^/0./  / 


0!(x,y)  EXP(-j47rf(y+x(jL>jjt)/c)  dxdy 


— /jW  —AW 

X  y 

where  (x)  =  rotation  rate 

R 


^  ^  e[  A(t^)A*(t2)]^  ffs(x,y)  EXP(-j47rf60j^x(t^-t2)/c)  dxdy 
E[|A(ti)|^]  Jj6(x,y)  dxdy 

Define  a  normalized  crossrange  scattering  function  S^(x)  as: 

s^(x)  ^ 

JJe(x,y)  dxdy 


Then  the  correlation  function  is  given  as: 

EXP(-J27r  l^x)  dx 

3  P(0)  -  jzTTl^  fxS  (x)  dx  -  27T^I/ 


( X )  dx 
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APPENDIX  A. 2  NONCOHERENT  INTEGRATION  OF  CORRELATED  GAUSSTAN(i.e.  RAYLE13' 
AMPLITUDE)  RETURNS. 


The  received  data  consists  of  M  complex  returns,  z  ,z, 

1  t. 


the  data  vector  Z  represent: 


The  components  7.^  are  samples  of  a  correlated  complex  gaussian  process.  The 
covariance  matrix  of  the  samples  is: 


J\.^  =  ^  [  Zj,  j 


where  the  operation  (  )'^  represents  conjugate  transpose. 


The  probability  density  function  of  receiving  a  specific  vector,  Z^^,, 
is  given  by  a  gaussian  density: 

P(Z„) .  (  (27r  )"IAJ  )-i  EXP(  -X  zihr^  Z,, ) 

It  is  desired  to  define  a  test  between  the  two  hypotheses  and  evaluat* 


the  performance  of  the  test. 

H  =  Hypothesis  0:  Data  is  noise  alone 


.  A..  =  '4  N. 


M  O 


=  Hypothesis  1:  Data  is  Rayleigh  signal  +  noise  so  that 
A.j^  -  /$  I  is  positive  definite. 

The  class  of  tests  to  be  evaluated  is  the  quadratic  forms  in  Z  : 


\ 
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Performance  of  this  test  can  be  found  if  the  probability  density  function 
of  q  can  be  derived.  This  can  be  accomplished  by  first  finding  its  charact¬ 
eristic  function. 


0„(CJ)  =  E[eJ^'^]  =  //.../EXP(y.  P(Zj,)  d^^d.^...dz,^ 

j-1  ij... ./eXP(-^ 
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viiere  the  limits  on  all  integrals  are  from  -  OO  to  +  OO . 


Let  be  a  unitary  transform  vAijch  diagonalizesAj^r^f,  » 
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where  is  the  m-th  eigenvalue  Qj^.  Note  that  I'yj  =  1  and 


I  here fore: 
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This  integration  can  be  carried  out  for  a  variety  of  cases.  Only  two 
cases  will  be  shown  here  assuming  simple  noncoherent  integration,  i.e. 

CASE  1:  Distinct  Eigenvalues,  Signal  +  Uttiite  Noise. 

3y  a  partial  fraction  expansion  one  obtains  that: 
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where  U(q)  is  a  unit  step  function. 
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The  probability  of  detection  is  then: 
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CASE  2:  M  Equal  Eigenvalues.  Noise  Alone. 
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Performance  of  this  test  can  be  found  if  the  probability  density  function 
of  q  can  be  derived.  This  can  be  accoaplished  by  first  finding  its  charact 
eristic  function. 

W  )  .  E  [  piz„) 

=  ((27r)"!A..Jr"  //.... /exp(-!s  z;,(A-;,-jWo„)z„  dp,...d. 
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vAiere  the  limits  on  all  integrals  are  from  -OOto  +  OO . 

Let  be  a  unitary  transform  v^ich  diagonalizesyVj^Qj,  * 

-  I  «nd  0„^  =  D„ 


where  is  the  m-th  eigenvalue  of.^^  Note  that  =  1  and 
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This  integration  can  be  carried  out  for  a  variety  of  cases.  Only  two 
cases  will  be  shown  here  assuming  simple  noncoherent  integration,  i.e. 


CASE  1:  Distinct  Eigenvalues.  Signal  +  White  Noise. 
3y  a  partial  fraction  expansion  one  obtains  that: 
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where  U(q)  is  a  unit  step  function. 
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The  probability  of  detection  is  then: 
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The  probability  of  false  alarm  is  then: 

TO  .TO 


fp^fq)  dq  =  (2/N^)''’  (l/r(M))  T  q”  ^ 
=  EXP(-27'y'N^)  \  (m!)~^  (2'y/N^)"' 


EXP(-2q/N  )  dq 


7^0 


m=0 


•/here  is  the  Gamma  Functionv  i.e.  I** fM)  =  (H-1)! 


dCJ 


73 


APPENDIX  A. 3  Dictr ibutioa  or  Monopuioo  statistics. 

A 

^  f  i 

'rile  suju  dn^  c  i I  -  tjx'OiiVi ;i  ^  c-t*e  d»ji'i‘  wlcibsu  eofipleic 

Kaussian  random  prooaasen,  Thi  jo-nb  '^rooabillty  densitji  firn^tion  of 
simultanco-a  aamplos  from  -ciie  t'Jo  pcjcesses  is  given  by: 
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A  change  of  variables  can  now  be  done.  Let; 
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Note  that  this  function  depends  only  on  -  (j)^  and  not  on  the  values  of 
the  angles  themselves.  Defining  ~  ~  integrating  oiit  one 

variable  yields: 

pCVi.V^.l//)  =  (27r)"^|y\.|"^Vj^V2EX?(-J^(/9^^V^+ ^22^2+2^12^2  ^ 

From  this  density  function,  the  density  of  the  monopulse  channel  statistics 
can  be  derived.  Defining  X  and  Y  as  the  real  and  imaginary  parts  respec¬ 
tively  of  the  monopulse  channel  then: 
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The  change  of  variables  from  )  to  (X,Y,V^)  has  Jacobian, 


J  =-  (V^/Vj). 


tlX.Y.V^  )  =  (27r(A|)"^V^  EXP(-1*V^(^^^  +  +  2^^^)) 

p(X,Y)  =  p(X,Y,V^)  dV, 
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p(X,Y)  = 


TT  (w‘'+m^+X%Y^-2mX)^ 


7r((X-H)^+Y^+W^)^ 


Def ini: ig  U  =  X  -  m; 


yr 


p(u,y)  = 


'TT  Mr-  .  ‘* 
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oo 
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.  O/'l 


dU  r=  0 


dY  =  0 


p(U|Y)  =  p(U,Y)/p(Y'  = 


2(W^+Y^)"*^ 
TT  (W^+Y^+U^)^ 
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oo 
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Data  will  be  thresholded  on  Y,  discarding  all  data  when  |Y(—  Y^.  The 
variance  of  tlie  monopulse  estimates  for  the  data  retained  is: 


i  1 


(Y^/W) 


APPENDIX  A. 4  High  SNR  Approximation  to 

2 

From  Section  3.2.1  the  width  parameter  W  is: 

/  =  -  ( P 

For  a  target  with  RCS  centroid  on  the  boresite  ft  =  0: 

5  n  SW  +  N  /2b^ 

=  Pvy  .  -L-, - 

/S(e>)  (J&  +  N^/2a 

With  hi^  SNR,  i.e.,  j SW  d!»:Si>  Ny2a^ 

„  S(«-)  de  N 

,  2  J  O 

W  —  +  ~  2  ^ 

I  S(«-)  dO  2b  /  S(©-)  c» 


Recognizing  that  the  first  term  is  ttie  noise  free  width  parameter, 

2 

W  ,  and  the  second  term  is  related  to  the  sum  channel  SNR: 

o  ’ 


SNR  = 


/  S(6-) 


Then  it  is  possible  to  approximate  W  as: 


+  (a^/b^)SNR"^ 
o 

Assume  a  quadratic  beam  shape,  f^(i)  =  1  -  1,17((J  -  » 

where  is  the  3-db  beamwidth  and  the  beams  cross  at  tiie  3-db  points. 
B  W 


a  =  f^(0)  =  0.707 


b  =  |f'(0)|  = 


=  W^  +  0.365  SNR"^ 
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APPENDIX  A. 5  Maximum  Likelihood  Estimation  of  Monopulse  Angle 
for  Raylei^  Scatterer. 

The  sum  and  difference  channels, X  and  ^ ,  are  correlated  gaussian 
random  variables  with  density  as  given  in  Appendix  A. 3  when  the  target 
is  a  Raylei^  scatterer.  Assume  that  M  pulses  have  been  transmitted  and 
the  data  consists  of  M  SEunples  of  each  channel,  ^  k=l,2,.,.M. 

Bach  data  sample  Is  complex  and  it  is  assumed  that  there  is  no  pulse-to- 
pulse  correlation.  The  joint  probability  density  function  of  receiving 
all  the  data  is; 


M 


p(L^Aj^;k=i.2,...M)  =  (27rP’IAh”Exp(-)^  ^  (Zk A' 


k=l 


■(£)’ 


A  =  4 


P 


When  there  is  no  noise,  the  matrix/ V  can  be  written: 

1 


A  =  ^ 


where  m  and  W  are  as  defined  in  Section  3.2.1  and  // 

S(^)  dfr  .  The  maximum  likelihood  estimate  of  m  is  the 
value  which  maximizes  the  probability  density  function.  That  is  one  must 
find  the  valu^  of  m  vdiich  satisfies  the  following  expression: 

<z:  -Vi; 

y  V^k> 
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=  0 


rm= 


III 
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i.e.  the  best  estimate  of  the  angle  given  many  independent  monopulne  angle 
estimates  is  tlie  RCS  weighted  sum  of  the  single  pulse  angle  estimates.  This 
also  demonstrates  that  for  a  single  pulse  of  data,  the  conventional  monopulse 
angle  estimate  is  optimim  under  the  assumptions  made. 

The  next  step  is  to  find  the  mean  and  variance  of  m: 
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Jsal _ 
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III. 


Since  are  gaussian  with  covariance  matrix 
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♦  These  equations  are  well  known  and  derived  from  the  orthogonality  prin. 
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Note  that  q  is  chi-square  order  2M. 
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p(q)  =  (  /i”  (M  -  1)1  )"^  q”'^  EXP(-q/  fJL^) 
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